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TECHNICAL  NOTE  3067 

ROLLING  EFFECTIVENESS  AND  AILERON  REVERSAL  OF 
RECTANGULAR  WINGS  AT  SUPERSONIC  SPEEDS 
By  John  M.  Hedgepeth  and  Robert  J.  Kell 

SUMMARY 


Linearized  supersonic  lifting-surface  theory  is  used  in  conjunction 
with  structural  influence  coefficients  to  formulate  a method  for  ana- 
lyzing the  aeroelastic  "behavior  in  roll  at  supersonic  speeds  of  a rec- 
tangular wing  mounted  on  a cylindrical  body.  Rolling  effectiveness  and 
aileron -reversal  speed  are  computed  by  using  a numerical  solution  which 
incorporates  matrices. 

Results  obtained  for  an  example  configuration  by  using  this  method 
are  compared  with  the  results  obtained  by  using  simplified  methods  of 
analysis.  For  the  particular  configuration  considered,  the  variation  of 
rolling  effectiveness  with  Mach  number  is  found  for  two  constant-altitude 
flights. 


INTRODUCTION 


In  the  past,  most  aeroelastic  calculations  have  been  based  on  the 
use  of  beam  theory  for  the  structural  analysis  and  strip  theory  for  the 
aerodynamic  analysis.  Hie  application  of  these  simplified  theories 
avoids  complications  which  result  from  using  more  refined  theories;  in 
addition,  the  simplified  theories  are  quite  adequate,  in  most  cases, 
when  applied  to  wings  of  high  aspect  ratio.  When  applied  to  wings  of 
low  aspect  ratio,  however,  these  simple  theories  may  become  inadequate; 
if  so,  more  refined  structural  and  aerodynamic  analyses  are  necessary. 

The  purpose  of  this  paper  is  to  describe  a method  for  predicting 
aeroelastic  effects  on  the  steady-state  roll  of  rectangular  wings  at 
supersonic  speeds  in  those  cases  for  which  beam  and  strip  theory  are 
inadequate  but  for  which  the  aerodynamic  effects  of  chordwise  deforma- 
tion may  be  neglected.  In  this  method  the  structural  distortions  caused 
by  arbitrary  loads  are  expressed  in  terms  of  structural  influence  func- 
tions. Hie  aerodynamic  loads  caused  by  arbitrary  angle -of -attack  distri- 
butions are  determined  by  superposing  basic  aerodynamic  loadings  resulting 
from  unit-step  angle -of -attack  distributions,  which  loadings  are  obtained 


2 


NACA  TN  3067 


herein  on  the  basis  of  linearized  aerodynamic  theory.  By  means  of  this 
superposition  procedure,  first  outlined  by  Frick  and  Chubb  in  refer- 
ence 1,  the  application  of  three-dimensional  lifting -surface  theory  is 
considerably  simplified. 

The  aircraft  configuration  considered  herein  consists  of  two  flex- 
ible rectangular  wings  with  trailing-edge  ailerons  of  constant  chord, 
diametrically  mounted  on  an  infinitely  long,  rigid,  cylindrical  fuselage. 

The  analysis  of  the  aeroelastic  rolling  behavior  is  separated  into 
various  parts.  The  analysis  of  the  structural  deformations  is  described, 
the  aerodynamic  loads  are  then  found,  and  the  two  parts  are  combined. 

A numerical  solution  of  the  resulting  equations  is  presented  in  matrix 
form.  Tables  of  aerodynamic  matrix  elements  usable  for  any  rectangular 
plan  form  are  included.  A particular  example  is  analyzed  and  the  results 
are  compared  with  those  obtained  by  simplified  methods.  Hie  structural 
analysis  of  the  example  configuration  is  included  in  appendix  A and  the 
details  of  the  aerodynamic  analysis  are  relegated  to  appendix  B. 


SYMBOLS 


D 

E 

G 

GM(y,n) 

L(y) 

M 


local  flexural  stiffness,  EtVl2(l  - \fi) 

Young's  modulus  of  elasticity 

shear  modulus  of  elasticity,  E/2(l  + p) 

structural  twist  influence  function  which  results  from  a 
unit  concentrated  load  at  the  wing  midchord 

structural  twist  influence  function  which  results  from  a 
unit  concentrated  torque 

aerodynamic  load  per  unit  span,  positive  upward 

free -stream  Mach  number 


M(y)  aerodynamic  moment,  per  unit  span,  about  the  midchord, 

positive  in  the  positive  twist  direction 


static  pressure  at  altitude 


Pq  standard  static  pressure  at  sea  level 

Q(y)  aerodynamic  moment,  per  unit  span,  about  the  elastic  axis, 

positive  in  the  positive  twist  direction 
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V 

a 

b 

t>a 

c 

ca 

e(y) 

l 

m 

P 

pb/2V 

q 

t 


x,y,z 


3 

3cz(y) 

3cm(y) 

3cq(y) 

6 

e(y) 

A 


free -stream  velocity 

ratio  between  fuselage  radius  and  exposed  wing  semispan 

total  wing  span,  2(a 1 + l) 

single  aileron  span 

wing  chord 

aileron  chord 

distance  measured  forward  from  the  midchord  to  the  elastic 
axis,  expressed  as  fraction  of  chord 

exposed  wing  semi span 

modified  aspect-ratio  parameter,  3z/c 

rolling  angular  velocity  (see  fig.  l) 

tangent  of  the  wing-tip  helix  angle 

dynamic  pressure 

thickness  of  wing  cross  section 

coordinate  system  (see  fig.  l) 


cotangent  of  the  Mach  angle. 


section  lift  coefficient. 


L(y) 

qc/p 


section  moment  coefficient  about  the  midchord. 


M(y) 


qc2/p 


section  moment  coefficient  about  the  elastic  axis, 

aileron  deflection  (see  fig.  l) 
angle  of  twist  of  wing  (see  fig.  1) 


Q(y) 

qc2/p 


structural  parameter,  — \j 2k(l  - p) 


j.  . o , , 

: ?* 


k 
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(i  Poisson's  ratio 

<r(x,y)  local  angle  of  attack  of  wing 

0 rolling  effectiveness,  (pb/2V)p/(pb/2V)p 

velocity  potential 

$(xi,yi)  nondimensionalized  velocity  potential 

Subscripts: 

F flexible  wing 

R rigid  wing 

a airloads  due  to  aileron  deflection 

r airloads  due  to  roll 

s airloads  due  to  structural  deformation 

rev  aileron  reversal 

p aerodynamic  coefficients  due  to  unit  pb/2V 

pQ  aerodynamic  coefficients  due  to  a unit  rate  of  roll  about 

the  x-axis 

a aerodynamic  coefficients  due  to  a unit  angle  of  attack  of 

the  entire  wing 

6 aerodynamic  coefficients  due  to  a unit  aileron  deflection 

0 indicial  aerodynamic  coefficients  due  to  a unit-step  angle - 

of-attack  distribution 

I, II, ...  VI  regions  on  surface  of  wing 


nondimensional  quantities  used  in  appendix  B 


Superscript: 


aerodynamic  coefficients  due  to  a unit-step  angle -of -attack 
distribution  on  one  wing  only 
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ANALYSIS 

Structural  Deformations 


In  this  report  the  treatment  of  the  structural  deformations  is 
based  on  the  assumption  that  there  is  no  chordwise  bending  of  the  wing. 
The  effect  of  the  structural  distortion  on  the  aerodynamic  loads  is 
then  determined  entirely  by  the  angle  of  twist  0( y),  the  contribution 
of  the  spanwise  bending  being  negligible.  (See  fig.  1.)  Consequently, 
only  the  determination  of  9(y)  will  be  included  in  this  analysis. 

If  the  section  lift  L(y)  and  the  section  moment  about  the  mid- 
chord M(y)  are  known,  the  angle  of  twist  can  be  obtained  from 


9(y)  = / GL(y,T])  l^)  dTl  + / Gm(y^)  M(q)  drj  (l) 

Jo  Jo 

In  this  equation,  Gj^(y , t])  and  G^(y,T))  are  influence  functions  which 

define  the  twist  at  y caused  by  the  application  of  a unit  concentrated 
load  at  the  midchord  and  a unit  concentrated  torque,  respectively,  at 
the  station  q.  For  many  structures  these  influence  functions  may  be 
obtained  analytically  as  is  done  in  appendix  A for  a uniform  flat-plate 
wing;  other  structures  may  be  handled  analytically  by  methods  such  as 
those  described  in  references  2,  3,  and  4,  for  instance.  For  some 
structures,  it  may  be  more  convenient  or  even  necessary  to  resort  to 
the  use  of  experimental  influence  coefficients.  The  analysis  proceeds 
hereinafter  on  the  assumption  that  the  influence  functions  are  known. 

In  general,  the  two  influence  functions  are  needed  in  order  to 
specify  completely  the  twist  of  the  wing.  For  many  rectangular  wings, 
however,  sufficient  accuracy  can  be  obtained  by  expressing  the  twist 
solely  in  terms  of  the  moment  about  some  "elastic  axis. " This  elastic 
axis  is  herein  defined  as  a line  along  which  loads  can  be  placed  without 
producing  significant  twist  anywhere. 

If  an  elastic  axis  does  exist,  it  is  no  longer  necessary  to  know 
the  influence  function  associated  with  load;  only  %(y,il)  need  be 
determined.  The  twist,  in  this  case,  is  given  by 


e(y) 


G^y,^  Q(n)  dq 


(2) 
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In  this  equation,  Q(y)  is  the  section  torque  about  the  elastic 
axis  and  can  be  expressed  in  terms  of  L(y)  and  M(y)  by  means  of  the 
simple  moment -transfer  relation 


Q(y)  = M(y)  - e(y)  c L(y)  (3) 

where  e(y)  is  the  distance  measured  forward  from  the  midchord  to  the 
elastic  axis,  expressed  as  a fraction  of  the  chord. 


Aerodynamic  Loads 

The  section  lift  L(y)  and  moment  M(y)  may  be  expressed  in  coef- 
ficient form  as 


Uy)  = y 3cz(y) 
M(y)  = ^ pc^Cy) 


00 


For  convenience,  the  products  of  (3  = Vw  - 1 and  the  coefficients  are 
considered  herein  rather  than  the  coefficients  themselves. 

It  is  assumed  that  the  aircraft  is  undergoing  a steady  roll  about 
the  axis  of  the  body  and  that  this  axis  is  in  line  with  the  direction 
of  flight.  Consequently,  the  resulting  loads  are  due  solely  to  struc- 
tural twist  of  the  wing,  the  rolling  velocity  itself,  and  the  deflection 
of  the  ailerons.  Since  linear  aerodynamic  theory  is  to  be  used,  the 
coefficients  resulting  from  this  steady  rolling  maneuver  can  be  expressed 
as 


eci  - + (Pcz)r  + (ecz)a 

Pcm  = (Pcm)a  + (Pcm)r  + (Pcm)a 


(5) 


where  s,  r,  and  a refer,  respectively,  to  structural  deformation, 
roll,  and.  aileron  deflection. 

In  the  determination  of  the  coefficients  on  the  right-hand  side  of 
equations  (5)  the  wing  and  body  must  be  treated  as  a unit.  In  view  of 
the  fact  that  in  this  problem  the  pressures  resulting  from  the  presence 
of  the  body  are  important  only  in  the  neighborhood  of  the  wing  root  and 
therefore  contribute  only  slightly  to  the  structural  distortions  and 
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rolling  moment,  a rather  simple  Idealization  Is  made  regarding  the 
body;  that  Is,  the  body  Is  replaced  by  a rigid  plate  Joining  the  two 
wing  roots  together.  The  width  of  this  plate  is  the  same  as  the  body 
diameter.  (See  fig.  2(a).)  While  it  is  not  argued  that  this  ideali- 
zation is  the  most  nearly  correct  one,  the  resulting  aeroelastic  model 
does  have  the  advantage  of  simplicity  and  allows  aerodynamic  interaction 
between  the  two  wings.  The  presence  of  the  cylindrical  body  is  further 
taken  into  account  by  neglecting  the  effects  upon  the  rolling  moment  of 
the  pressures  acting  on  the  rigid  plate,  since  the  pressures  acting  on 
the  actual  cylindrical  body  would  produce  no  rolling  moment. 


Although  the  details  of  the  calculations  of  the  aerodynamic  loads 
have  been  relegated  to  appendix  B,  a short  discussion  of  the  loads  due 
to  each  of  the  three  causes  - structural  deformation,  roll,  and  aileron 
deflection  - is  included  in  the  subsequent  portions  of  this  section. 


Loads  due  to  structural  deformation. - As  a consequence  of  the  pre- 
viously assumed  linearity  of  chordwise  deformation,  the  local  angle  of 
attack  of  the  wing  is  completely  defined  by  the  twist  of  the  wing.  Since 
the  twist  is  not  defined  explicitly,  the  deformations  being  dependent  on 
loads  which  are,  in  turn,  dependent  on  the  twist  itself,  it  is  necessary 
to  be  able  to  perform  the  rather  difficult  task  of  calculating  the  aero- 
dynamic loads  caused  by  an  arbitrary  angle -of -attack  distribution.  This 
task  is  considerably  simplified  for  the  rolling  problem  by  superposing 
loads  caused  by  an  antisymmetrical  unit-step  angle -of -attack  distribution 
obtained  by  imposing  a positive  unit  angle  of  attack  outboard  of  any 
spanwise  station  q on  the  right  wing  and  a negative  unit  angle  of 
attack  over  the  corresponding  portion  of  the  left  wing  as  shown  in  fig- 
ure 2(b).  These  basic  load  distributions,  which  have  the  nature  of  aero- 
dynamic influence  functions,  are  hereinafter  called  "indicial"  loads  for 
the  sake  of  brevity.  Superposition  of  the  indicial  section  coefficients 
of  lift  and  twisting  moment,  designated  pcjQ(y,q)  and  Pc^^q), 


respectively,  which  result  from  such  basic  angle -of -attack  distributions, 
yields  the  section  coefficients  [pcj(y)JB  an<*  [Pcm(y)]6  t*ue  8^'ruc_ 


tural  deformation.  The  required  superposition  integrals  are  given  by 


°f0l  !? d" 
^h,n)  | 4" 


(6) 


where  the  twist  at  the  root  is  assumed  to  be  zero. 


■*v 
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The  indicial  section  coefficients  are  dependent  not  only  on  the 
spanvise  coordinate  y and  the  position  of  the  step  ti  but  also  on 
two  additional  parameters.  These  parameters  are  the  modified  aspect- 
ratio  parameter  0 z/c,  which  usually  appears  in  theoretical  supersonic 
aerodynamic  calculations,  and  the  nondimensional  body  radius  a.  The 
dependency  on  the  body  radius  is  undesirable  since  it  restricts  the 
application  of  these  indicial  section  coefficients  to  a particular  value 
of  a.  Fortunately  the  dependency  can  be  eliminated  by  making  use  of 
the  fact  that  the  loads  due  to  the  antisymmetrical  unit-step  angle-of- 
attack  distribution  can  be  separated  into  two  parts:  the  first  is  the 

loads  due  to  a unit-step  angle  of  attack  on  the  right  wing  only  (see 
fig.  2(c));  the  second  is  the  loads  due  to  a negative  unit-step  angle 
of  attack  on  the  left  wing  only.  It  can  be  seen  that  the  second  part 
is  merely  the  negative  mirror  image  of  the  first  and,  consequently, 
only  the  first  case  need  be  considered  in  detail.  The  total  indicial 
section  coefficients  can  be  written  in  terms  of  these  partial  indicial 
section  coefficients,  designated  0cZq'  and  00^',  as  follows: 


3cz0(y,Tl)  = 0cz0'(y,n) 
3cmQ(y>Tl)  = 0Cme,(y,n) 


0cz  '(-2al-y,T)) 


3cmg ' (-2al-y,q) 


(7) 


The  first  term  on  the  right-hand  side  in  equations  (7)  gives  the  contri- 
bution of  the  right-hand  step  angle  of  attack;  the  second  term  gives  the 
contribution  of  the  negative  left-hand  step. 


The  advantage  of  the  foregoing  separation  is  that  3cZq'  and 


are  independent  of  the  body  radius  and  are  functions  of  only  the  modi- 
fied aspect-ratio  parameter,  provided  that  this  parameter  is  greater 

than  - — that  is,  when  the  Mach  number  is  great  enough  so  that 

there  is  no  point  on  either  wing  that  is  influenced  by  both  wing  tips 

simultaneously.  Thus,  if  the  restriction  ~ > j '~2&  *s  imPosed>  a 

single  parameter  remains  - 0z/c  itself  - and  it  is  feasible  to  compute 
tables  applicable  for  any  body  radius.  Numerical  values  of  0 czq' 


and  0cfflQ ' 


for  - -fL  £ £ % 1 and  0 % ^ % 1 
pl  I X 


in  steps  of  0.1  are  pre- 


sented in  table  I for  several  values  of  0Z/c.  These  numerical  values 
were  obtained  from  expressions  derived  in  appendix  B. 
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Loads  due  to  roll. - 
from  roll  can  be  written 


Hie  section  loading  coefficients  resulting 
in  the  form 


(*l>r  - *lp  S 

(0cm)r  = Pcmp  fv 


(8) 


The  rolling  derivatives  pc^  and  P°mp  are  &^ven  in  appendix  B. 

Here  again  the  derivatives  are  dependent  on  both  pi/c  and  a.  Inspec- 
tion of  the  resulting  expressions  (eqs.  (B13)  and  (Bl^)),  however, 
reveals  that  the  dependence  on  a is  quite  simple;  the  expressions 
for  pc j and  pc^  can  each  be  separated  into  two  parts  - one  with 


the  coefficient  — - — , the  other  with  - both  of  which  parts 

1 + a*  1 + a 

are  independent  of  a.  It  is  possible  to  show  that  in  each  case  the 
first  part  is  merely  the  derivative  which  results  from  a unit  angle  of 
attack  of  the  entire  wing  and  the  second  part  is  the  rolling  derivative 
which  would  result  if  the  wing  were  rolling  about  the  x-axis.  Hiere- 
fore,  these  expressions  become 


0C7  = 


0C. 


"P 


1 + a PCla  + 1 + a PClPc 


1 + a PCnk  + 1 + a PCmP„ 


(9) 


In  order  for  this  simplification  to  be  correct,  it  is  again  neces- 
sary to  impose  the  restriction  — > — i . 

c 1 + 2a 

The  quantities  (3c,  , (3c_  , pc,  , and  Pc_  have  been  calcu- 

a “'a  fcP0  “P0 

lated  for  values  of  0 s j = 1 in  steps  of  0.1  and  are  presented  in 

ir 

table  II  for  several  values  of  pz/c. 


■ -ir-  :*>  “I  . 


?fr  ' 1 ’ 
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Loads  due  to  aileron  deflection. - The  section  coefficients  of  lift 
and  tvistlng  moment , 


and 


(0C z)a  = Pcl66 

(pcn^  = pc^b 


(10) 


respectively,  are  found  for  constant-chord  trailing-edge  ailerons  of 
arbitrary  length.  All  gaps  between  the  wing  and  aileron  are  considered 
to  be  sealed.  The  calculation  of  the  aileron  derivatives  PCj  and 


PCj^  is  very  similar  to  that  of  the  indicial  structural  loads  and  is 

included  in  appendix  B.  For  most  reasonable  aileron  configurations,  if 
the  aircraft  is  flying  at  supersonic  speeds  sufficiently  high  to  satisfy 

/ 8 2 1 

the  Mach  number  limitation  previously  imposed  [that  is,  -E—  > 


\ ' c 1 + 2a 

an  aileron  deflection  on  the  left  wing  produces  no  loads  on  the  right 
wing.  For  this  reason,  only  the  case  is  considered  wherein  the  loads 
on  one  wing  are  independent  of  the  aileron  deflection  on  the  other.  In 
this  case,  the  limitation  on  the  modified  aspect-ratio  parameter  pz/c 
for  the  analysis  of  aileron  loadings  in  appendix  B is  that  pz/c  must 

c / c 

be  greater  than  — — . Numerical  values  are  given  in  table  III 


1 + 2a  - 


ba* 

Z 


for  0 % y % 1 in  intervals  of  0.1  for  several  values  of  pz/c;  the 
computations  have  been  made  for  ba/z  = 1.0  and  ca/c  = 0.2. 


Significance  of  Mach  number  limitations. - Although  the  restrictions 
that  have  been  placed  upon  the  modified  aspect-ratio  parameter  limit  the 
utility  of  the  aerodynamic  results  contained  herein,  these  restrictions, 
in  reality,  are  not  serious.  This  fact  is  substantiated  by  considering 
a typical  configuration  such  as  that  used  for  the  example  contained  in  a 
subsequent  section.  For  this  wing,  which  has  an  aspect  ratio  of  3*6, 
the  pertinent  parameters  are  z/c  = 1.5,  a = 0.2,  ca/c  = 0.2,  and 

ba/z  = 1.  The  restriction  imposed  upon  the  expressions  for  the  loads 

due  to  structural  deformation  and  roll,  — > — - — , and  that  for  the 

' c l + 2a' 

loads  caused  by  aileron  deflection,  become  M > 1.108 

1 + 2a  - -^ 

and  M > 1.05^,  respectively.  Since  the  validity  of  linear  aerodynamic 


>*$*;*.:  i"? 
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theory  Is  questionable  near  a Mach  number  of  1,  these  limitations  are 
of  little  consequence. 


Aeroelastic  Solution 


Structural  and  rolling-moment  equations.-  In  the  solution  of  the 
aeroelastic  rolling  problem,  not  only  must  structural  equilibrium 
(eq.  (l))  be  satisfied  but  also  the  equilibrium  of  moments  about  the 
rolling  axis.  If  the  loads  in  equation  (l)  are  replaced  by  the  loads 
arising  from  the  various  causes  as  derived  in  the  preceding  section, 
the  equation  specifying  structural  equilibrium  becomes 


- t/0 


r ^ d9 

GL(y,n)  J 3cl0(n,£)  ££  + 


OM(y.n)/ol  ecmjOi.O  f§  <«J 


dr)  + 


qc  pb  r ^ 

— 2vJ0  pcipM  + c dTi + 


Rolling-moment  equilibrium  is  attained  by  setting  the  total  rolling 
moment  equal  to  zero.  This  condition  can  be  written  as 


f (aZ  + tj)  L(tj)  drj  - 0 

J0 


Again,  the  results  of  the  preceding  section  can  be  used  to  give 


^ Jq  + ^ Jc\  + 57  *ipW  + & 3cz&(ti) 


dq  = 0 


(12) 


* 

2<?T 

S,  *2XULSLj&- 


a 


■Ip 
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Alternate  structural  equation. - Equations  (11)  and  (12)  completely 
express  the  necessary  conditions  for  this  aeroelastic  problem.  It  is 
to  be  noted  that  the  quantities  5 and  qc/p  are,  in  general,  known, 
the  structural  twist  0(y)  is  unknown,  and  the  aeroelastic  rolling 
rate  pb/2V  is  the  quantity  that  is  desired.  The  simultaneous  solu- 
tion of  these  equations  in  closed  form  is,  in  actuality,  extremely 
difficult  if  not  impossible  to  obtain.  Some  sort  of  numerical  solution 
is  therefore  inevitable.  One  type  of  numerical  solution,  which  utilizes 
a collocation  technique  in  the  solution  of  equations  (ll)  and  (12),  is 
derived  in  a subsequent  section. 

The  numerical  methods  with  which  the  indicated  integrations  in 
equations  (ll)  and  (12)  are  performed  can  be  expected  to  be  accurate 
when  applied  to  reasonably  well-behaved  functions.  In  addition  to 
numerical  integration,  however,  numerical  differentiation  must  be  used 
in  order  to  express  the  rate  of  twist  de/dy,  upon  which  the  loads  due 
to  structural  deformation  depend,  in  terms  of  the  twist  9.  As  is  well 
known,  the  process  of  numerical  differentiation  is  not  as  accurate  as 
that  of  numerical  integration  and  therefore  should  be  avoided  if  possible. 
The  numerical  differentiation  can  be  eliminated  by  differentiating  equa- 
tion (ll)  with  respect  to  y,  which  operation  yields 


d9  _ qc  / 

dy  P Jo 


dGL(y,n)  pi  de  .. 

-V— J0  d?a;  + 


^(y,n)  rl  de 

-3T- j0  CW,,’C)«d{; 


drj  + 


In  this  manner,  equilibrium  equation  (ll)  is  replaced  by  equation  (13) 
and  the  twist  9 no  longer  appears  in  the  problem.  Hie  rate  of  twist 
becomes  the  unknown  and  can  be  used  directly  in  the  numerical  solution. 


The  use  of  equations  (12)  and  (13)  is  particularly  desirable  when 

aGL(y,Ti)  and  dGM(y>n) 

dy  dy 


the  functions 


-,  which  are  merely  the  rates  of 
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twist  due  to  concentrated  loads  and  torques,  respectively,  can  be  found 
accurately.  Hie  example  considered  in  a subsequent  section  is  one  for 
which  it  was  possible  to  obtain  these  functions  analytically.  When 
complicated  structures  are  to  be  dealt  with,  however,  recourse  must 
usually  be  made  to  experiment  or  to  approximate  theories  in  obtaining 
the  structural  properties.  For  these  structures,  probably  only  discrete 
influence  coefficients  for  twist  can  be  found}  the  determination  of  the 
rate-of -twist  influence  functions  would  be  subject  to  at  least  the  same 
inherent  inaccuracies  as  the  numerical  differentiation  discussed  pre- 
viously. In  such  cases,  therefore,  the  use  of  equation  (13)  instead  of 
equation  (ll)  would  not  be  advantageous. 

Aileron  reversal.-  If  the  rolling -moment  equation  (12)  is  solved 
for  5,  the  result  is  substituted  into  equation  (ll)  or  equation  (13), 
and  pb/2V  is  set  equal  to  zero,  the  following  equations  result: 

For  equation  (ll). 
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For  equation  (13) , 


<19  (y)  = <^revc 
dy  p 


r'fioLto*)  rl  „ / dp  J(. 

Jo  sr-Jo  « dC 


aoM(y.’i)  r»  d9 

ST~ J0  pc"e(,,t)  if  E 


dn  - 


/* 


aGL(y,Ti)  « / x a , , 

—5—  ec,8(,)  + 0 — 5^—  ^(1) 

(al  + n)  Pols(l)  dT| 


dq 


(aZ 


+ n)  /'  | '‘"I 


(15) 


Either  of  these  homogeneous  Integral  equations  expresses  the  condition 
for  aileron  reversal,  where  the  dynamic  pressure  at  reversal  qrev 

appears  as  an  eigenvalue. 

Rolling  effectiveness.-  Hie  effect  of  the  elasticity  of  the  wing 
on  the  rolling  behavior  may  be  determined  by  examining  the  ratio  between 
the  rate  of  roll  of  the  flexible  wing  and  the  rate  of  roll  which  would 
occur  if  the  wing  were  rigid.  Hiis  ratio,  the  "rolling  effectiveness," 
is  designated  as  f>,  where 

, _ (pb/2V)F 
P " (pb/2V)R 

In  this  equation  (pb/2V)F  is  the  rate  of  roll  determined  directly 
from  the  aeroelastic  equations.  The  quantity  (pb/2V)^,  the  rate  of  roll 
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for  the  rigid  wing,  is  found  by  setting  the  twist  per  unit  length  equal 
to  zero  in  equation  (12): 


f (a*  + y)  pc,  (y)  dy 

^0 o 

(a l + y)  3clp(y)  dy 


(16) 


Simplifications  for  wings  with  an  elastic  axis . - If  the  wing  is 
constructed  in  such  a manner  that  an  elastic  axis  does  exist  as  previ- 
ously described,  the  expression  for  structural  equilibrium  is  given  by 
equation  (2)  instead  of  equation  (l).  It  is  to  be  remembered  that,  in 
this  case,  the  twists  are  dependent  only  on  Q(y),  the  torque  about  the 
elastic  axis,  and  not  on  the  load  L(y).  The  aeroelastic  structural- 
equilibrium  equation  which  results  from  using  equation  (2)  instead  of 
equation  (l)  can  be  obtained  by  merely  deleting  the  terms  in  equa- 
tion (11)  that  involve  (^(y^)  and  replacing  the  section  moment  coef- 


ficients about  the  midchord  by  the  section  moment  coefficients  about  the 


elastic  axis.  Thus,  the  derivatives  pcm  , 

0 


Pcmp>  60(1 


0c; 


mg 


should  be 


replaced  by  pc 


0cv 


and  0c^,  respectively,  where 


pc 


*0 


3cm„  - e pc, 


PCqp  = Pcmp  - e pcZp 
0C% ' ' e 


(17) 


In  these  equations,  as  before,  e is  the  distance  measured  forward  from 
the  midchord  to  the  elastic  axis,  expressed  as  a fraction  of  the  chord. 


The  equations  analogous  to  equations  (lj),  (l4) , and  (15)  are 
obtained  in  exactly  the  same  manner. 


I 


NUMERICAL  METHOD 


Since  an  exact  solution  of  the  aeroelastic  equations  is  not  feasible, 
even  for  the  simplest  configuration,  this  section  is  devoted  to  the 


A 

H, 


! - ff 


O 
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presentation  of  a numerical  method  of  solution,  lhe  method  is  based  on 
satisfying  the  equations  at  a number  of  discrete  spanwise  stations. 
Matrix  notation  is  used  as  an  aid  in  organizing  the  numerical  procedure. 
A step-by-step  outline  of  the  procedure  is  included  at  the  end  of  this 
section  in  order  to  aid  the  reader. 

The  accuracy  of  the  numerical  method  is  directly  dependent  upon  the 
number  of  stations  used.  Experience  has  shown  that  for  this  particular 
problem  11  stations,  defined  by  the  end  points  of  10  equal  spanwise 
intervals,  are  sufficient.  Hie  derivation  that  follows  is  therefore 
based  on  10  equal  intervals,  the  extension  to  other  numbers  of  intervals 
being  evident.  Simpson’s  rule  is  used  to  perform  the  integrations  and 
parabolic  difference  equivalents  are  used  to  replace  any  necessary 
derivatives. 


Matrix  Operations 

Two  distinct  types  of  integration  appear  in  the  aeroelastic  equa- 
tion. The  first  is  of  the  form 


and  the  second  is  of  the  form 


h(y,n)  g(q)  dq 


Hie  integrands  in  these  equations  are  evaluated  at  each  of  the  spanwise 
stations  t)q,  t^,  tj2,  . . . 7^,  . . . t)10,  where,  for  equal  intervals  of 

width  e,  = ie.  In  a similar  manner  this  subscript  notation  is  used 

to  denote  evaluation  of  the  integrands  at  each  spanwise  station.  Thus, 
fi  = f(’li)  = f(ie),  gi  = g^i)  « gUO,  and  hy  = hfy^Tij)  = h(ie,je). 
For  10  intervals,  e = l/ 10  and  Simpson's  rule  becomes,  for  each  of  the 
integrals,  respectively, 


J f(n)  s(n)  dTj  = ^(foSo  + 4flSl  + 2f262  + 

^f3®5  + • • • + ^989  + flo8lo) 


& 


■i 

1 ,1  - 


■it- 
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: \ ' 


and  (for  any  single  y^) 

f0  h(yi^)  g(n)  dr]  = 5o(hi0g0  + ^hilgl  + ^±2*2  + 

^1383  + • • • + ^1989  + hi  ,ioeio) 

In  matrix  form  these  integrals  can  be  written  as,  respectively, 

/0‘  f(n)  s(i)  dr,  = [fJOQkl 

and  (for  all  y^) 


The  row  matrix  is  given  by 

LfiJ  = [jo  fl  f2  f3  ***  f9  fioJ 

and  a typical  column  matrix  by 

81 

82 

s3 

kl  - • 

I I • 

*9 

sio 


18 


NACA  IN  5067 


where  the  subscripts  denote  the  position  of  each  element.  The  square 
matrix  is  given  by 


**00 

h01 

*02 

**0,10 

h10 

*11 

h12 

•••  hl,10 

0 

CVI 

A 

h2i 

h22 

•••  *2,10 

• 

• 

• 

• • • • 

h10,0  h10,l  h10,2  *•*  h10,10 


where  the  first  subscript  designates  the  row  and  the  second  subscript 
designates  the  column  in  which  an  element  appears.  This  subscript  nota- 
tion is  used  throughout  this  section.  The  integrating  matrix  is  given 
by 


(18) 


All  derivatives  which  appear  in  the  aeroelastic  equation  are  of 
first  order;  difference  equivalents  based  on  passing  a parabola  through 
three  adjacent  points  are  used  herein  to  approximate  the  derivatives. 
For  the  points  i * 1,  2,  . . .9  the  standard  difference  equivalent 
derived  by  finding  the  slope  at  the  center  of  the  three  ordinates  is 
used: 


fl+l  ~ fi-l 
2c 


T 


. . 


.-)■  x. 
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For  the  end  points  i ■ 0 and  i ■ 10,  the  slope  at  the  exterior  of 
the  three  ordinates  must  be  used! 


/df  \ ^ -3fp  + 4fj  - t2 
W/0  * 2e 


/df  \ _ f8  ~ ^f9  + ^f10 

W/10  * 2e 


In  matrix  form  with  e = l/ 10,  these  expressions  become 


where  the  differentiating  matrix 


-3  4-1 
-10  1 
-10  1 
-10  1 
-10  1 
-10  1 
-10  1 
-10  1 
-10  1 
-10  1 
1-4  3 


(19) 


Aeroelastic  Matrix  Equation 

The  aeroelastic  equations  may  now  be  written  in  matrix  form  by 
using  the  results  of  the  foregoing  numerical  analysis  of  integrating 
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and  differentiating  processes.  Equations  (11)  and  (12)  become, 
respectively. 
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If  these  equations  are  divided  through  by  5 and  the  resulting 
matrix  equations  are  combined,  the  following  partitioned  matrix  equa- 
tion results: 


m-i4 

0 I 

H 

1 0k/&| 

S 

01 1 

lAJ  ! 

E 

§/s 

P 

(22) 


where 


is  the  unit  matrix. 


The  submatrices  are  given  by 
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where 


i,  J,  k - 0,  1,  2,  . . . 10 


For  a given  configuration  flying  at  a given  speed  and  altitude, 
all  the  quantities  in  equation  (22)  are  known  except  the  twist  9 at 

each  station  and  the  rate  of  roll  solution  of  this  matrix  equa- 

tion yields  these  quantities. 

The  rigid  rate  of  roll,  as  found  by  expressing  the  integrals  in 
equation  (16)  in  numerical  form,  is  merely 


(25) 


If  the  rate-of -twist  influence  coefficients  are  known,  it  is  prefer- 
able to  employ  equation  (13)  instead  of  equation  (ll).  In  this  case,  the 
partitioned  matrix  equation  analogous  to  equation  (22)  is 


1 H&K 

M .[ 

L J L 

Aik 

W| 

1 H*  1 

1 

/d9\ 

W/k 

6 

= 

Ci 

1 

1* 
1 

J" 

E 

&I& 

2V/ 

F 

where  the  submatrices  are  given  by 


(24) 


?' 
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and  the  expressions  for  the  scalars  E and  F remain  the  same  as  those 
given  immediately  after  equation  (22). 

By  applying  the  same  numerical  processes,  the  aileron-reversal  equa- 
tion (14)  (or  eq.  (15))  can  be  put  in  matrix  form.  Equation  (14),  in 
which  0 is  the  variable,  becomes 


<lrevc 


[%k]  + |ci|  lOJ 


ek 


(25) 


The  submatrices  are  defined  immediately  after  equation  (22). 


H fit ' * 


/ 
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Similarly,  equation  (15),  where  dQ/dy  is  the  variable,  becomes 


The  submatrices  are  defined  immediately  after  equation  (24). 

A solution  to  equation  (25)  (or  eq.  (26))  is  easily  found  by  using 
matrix  iteration.  The  process  converges  to  give  the  lowest  eigenvalue 
from  which  the  dynamic  pressure  at  aileron  reversal  may  be  obtained. 


Reduced  Matrix  Equations 


Up  to  this  point,  the  numerical  analysis  has  been  based  on  the  use 
of  11  spanwise  stations,  a number  that  was  deemed  necessary  in  order  to 
obtain  the  desired  accuracy  because  of  the  ill-behaved  nature  of  some 
of  the  aerodynamic  loading  functions.  In  most  cases,  however,  the  twist 
(or  rate  of  twist)  is  well  behaved  and,  therefore,  requires  fewer  sta- 
tions for  adequate  specification.  Utilization  of  this  fact  allows  a 
considerable  saving  in  the  amount  of  work  necessary  to  solve  the  matrix 
equations  because  of  the  fewer  degrees  of  freedom  involved. 


If  either  the  twist  or  the  rate  of  twist  is  specified  at  the  even- 
numbered  stations,  an  interpolation  procedure  can  be  used  to  determine 
the  values  at  all  11  stations.  Ihe  particular  type  of  interpolation 
used  herein  is  obtained  by  passing  a fifth-degree  polynomial  through 
the  even-numbered  stations  and  then  evaluating  this  polynomial  at  the 
odd-numbered  stations.  Bais  interpolating  procedure  can  be  written  in 
matrix  form  as 


where  k = 0,  1, 
polating  matrix 


. 10  and  l = 0,  2,  4,  . 
is  given  by 


10  and  the  inter- 
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The  twist  (or  rate  of  twist)  now  needs  to  be  known  only  at  the 
even-numbered  stations,  the  quantities  at  the  odd-numbered  stations 
being  obtained  by  interpolation.  Therefore,  only  the  even  rows  of  the 
influence-function  matrices  are  necessary;  these  matrices  become  rec- 
tangular with  11  columns  and  6 rows.  When  these  simplifications  are 
used  to  write  the  matrix  equations  for  equations  (ll)  and  (12)  and  the 
resulting  equations  are  combined,  the  following  partitioned  matrix  equa- 
tion results: 
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where  the  submatrices  are  given  by 


Note  that  the  subscripts  j and  k assume  all  values,  both  even  and 
odd,  and  the  subscripts  m and  l assume  only  the  even  values. 
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Hie  analogous  partitioned  matrix  equation  to  be  used  when  the  rate- 
of -twist  influence  functions  are  known  is 


~k\z  1 ■ H 

V 

/de\ 

6 

- 

Cm* 

(29) 

L v J 

E 

2 VI 

F 

where  the  submatrices  are  given  by 
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The  scalars  E and  F are  given  after  equation  (28).  Here  again  the 
subscripts  j and  k take  on  all  values,  whereas  m and  2 assume 
only  even  values. 

The  aileron-reversal  equations  may  also  be  expressed  in  this 
reduced-order  form.  Thus,  equation  (14),  in  matrix  form,  is 


The  submatrices  are  defined  after  equation  (28).  Similarly,  equation  (15) 
becomes 


The  submatrices  in  this  case  are  given  immediately  after  equation  (29). 
Again,  standard  matrix  iteration  procedure  may  be  used  to  solve  equa- 
tion (30)  or  (31),  the  amount  of  work  being  approximately  one-fourth  as 
much  as  is  involved  in  the  iteration  of  equation  (25)  or  (26). 


Simplifications  for  Wings  With  an  Elastic  Axis 


Two  operations  are  required  to  modify  the  aeroelastic  matrix  equa- 
tions when  an  elastic  axis  exists.  First,  all  terms  involving  the 


are  deleted.  Second,  the  section-moment- 


r 1 

8gl 

1 

1? 
1 

or 

_dy 

matrices 
coefficient  matrices 


0c 


*0 


PcqT 


, and 


% 


|3c. 


"P 


, and 


0c, 


®6 


are  replaced  by 


in  accordance  with  equation  (17). 
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Computational  Procedure 


As  has  been  previously  mentioned,  the  task  of  calculating  the 
rolling  effectiveness  for  the  given  configuration  flying  at  a given 
altitude  and  speed  is  straightforward ; only  the  solution  of  a set  of 
simultaneous  equations  is  necessary.  If  a complete  knowledge  of  the 
aeroelastic  behavior  over  large  ranges  of  speed  and  altitude  is  desired, 
however,  the  manner  in  which  the  various  parameters  are  involved  should 
be  taken  into  account  in  determining  the  sequence  of  calculations.  It 
should  be  noted  that  the  Mach  number  enters  the  problem  in  a complicated 
manner,  whereas  the  dependence  on  altitude  is  rather  simple;  that  is, 
the  Mach  number  affects  the  aerodynamic  matrices  and  the  altitude  affects 
only  the  dynamic  pressure  q.  It  is  therefore  obvious  that  the  most 
economical  way  to  perform  the  computations  is  to  calculate  the  variation 
of  rolling  effectiveness  with  dynamic  pressure  for  several  constant  values 
of  Mach  number.  In  addition,  since  the  rolling  effectiveness  is  of  little 
interest  when  the  controls  are  reversed,  the  range  of  dynamic  pressures 
should  be  restricted  to  values  less  than  qrev.  The  dynamic  pressure  at 

reversal  should,  therefore,  be  determined  for  each  value  of  Mach  number 
before  proceeding  with  the  calculation  of  rolling  effectiveness. 


An  outline  of  the  steps  required  in  the  determination  of  the  aero- 
elastic effect  on  roll  for  a range  of  altitude  and  Mach  number  is 
included  herein.  For  simplicity,  only  one  of  the  several  numerical 
approaches  derived  in  this  section  - that  is,  the  one  wherein  the  twist 
influence  functions  G^  and  Gy  and  the  interpolation  procedure  are 

utilized  - is  illustrated;  the  others  follow  the  same  outline,  differing 
only  in  detail. 

(l)  Evaluate,  either  analytically  or  experimentally,  the  twist 
influence  functions  Gj/y,^  and  Gy(y,Tj)  at  stations  y/z  = 0,  0.2, 

0.4,  . . . 1.0  due  to  loads  and  torques  applied  at  stations  tj/z  = 0, 
0.1,  0.2,  . . . 1.0.  From  these  values,  form  the  ll-by-6  matrices 

[(=!■)  J 


(2)  For  a given  Mach  number  (one  should  be  chosen  which  results  in 
a value  of  pz/c  which  appears  in  the  tables)  determine  the  indicial 

section  lift  and  moment  coefficients  pcz  and  pc™  . Ihese  coeffi- 

9 v) 

cients,  which  are  found  by  applying  equation  (7)  to  the  values  of  pc  ’ 
and  pc'  in  table  I,  can  then  be  used  to  form  the  square  matrices 


:4-  -v.  ft- 


7 
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(3)  In  a similar  manner,  form  the  column  matrices 


(%)j 


and 


, making  use  of  equation  (9)  in  conjunction  with  the  values  of 


pc 


ZPr 


and 


PoIa 


in  table  II. 


(4)  Determine  the  section  loading  coefficients  due  to  a unit  aileron 
deflection.  Use  these  coefficients  to  construct  the  column  matrices 


(pcj&) 


and 


(These  coefficients  are  listed  for  one 


aileron  configuration  in  table  III.  Equations  for  the  coefficients  for 
a rather  general  aileron  configuration  are  included  in  appendix  B. ) 

(5)  Compute  the  matrices  )v|»  | cm*  | * and  [_  d2*J  and 

the  scalars  E and  F defined  immediately  after  equation  (28) . The 
integrating  and  differentiating  matrices  and  QdJ  in  these  defi- 

nitions are  given  by  equations  (18)  and  (19) i the  interpolating  matrix 
JjTkz]  Is  given  by  equation  (27) ; the  row  matrix  j^aZ  + yjj  is  made  up 
of  the  moment  arms  about  the  rolling  axis. 


(6)  Form  the  matrix 


M + 1 


Cm* 


(7)  Obtain  qrev  by  iterating  this  matrix  (see  eq.  (30)). 

(8)  For  each  of  several  values  of  q between  zero  and  qrev,  form 
the  matrix  equation  (28).  Solution  of  this  equation  yields 

for  each  value  of  q.  Obtain  the  rolling  effectiveness  0 by  dividing 

(l/s)r  * (£A)„  ((#)„  - 5} 

(9)  Repeat  steps  (2)  to  (8)  for  several  other  values  of  Mach  number. 


SAMPLE  APPLICATION 


In  this  section,  aileron  effectiveness  and  reversal  speed  are  found 
for  a specific  aircraft  by  the  approach  set  forth  in  this  paper,  which  is 
termed  "lifting  surface  theory,"  and  by  two  simplified  methods.  Each  half 
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of  the  exposed  wing  consists  of  a uniformly  thick  aluminum  liquate  with 
an  aspect  ratio  Z/c  of  3/2  and  a thickness  ratio  t/c  old.  .-02.  The 
full-span  aileron  is  formed  by  bending  the  plate  along  the  chord 

line}  therefore  ba/z  = 1.0  and  ca/c  = 0.2.  The  wings  ars mamounted 
rigidly  on  the  body,  which  has  a radius  of  one-fifth  the  erased  wing 
semispan;  therefore  a = 0.2.  These  ratios  are  sufficient  to  define 
completely  the  configuration,  since,  as  is  seen  later,  the  resasuJlts  are 
independent  of  absolute  dimensions. 

Although  the  plate  is  considered  to  be  bent  in  order  tif-Torm  the 
aileron,  the  plate  is  assumed  to  behave  structurally  as  if  nov  bend  had 
been  made.  On  the  basis  of  this  assumption  the  necessary  stnu.j-xc'tural 
influence  coefficients  are  determined  in  appendix  A by  meant  ocof  a sim- 
plified flat-plate  theory.  In  this  case,  as  could  be  expected^ EL,  not  only 
does  an  elastic  axis  exist  but  also  this  elastic  axis  coincllenas  with  the 
midchord. 


Calculation  by  Lifting-Surface  Theory 


Since  the  influence  functions  are  obtained  analytically  It  "tkie  rate- 
of -twist  influence  function  rather  than  the  twist  influence  [anj_nction  has 
been  determined  in  order  that  no  numerical  differentiation  leri  necessary. 
In  addition,  the  wing  structure  being  free  of  discontinuities,  "the  inter- 
polation procedure  developed  in  the  preceding  section  may  bent;  uccessfully 
employed.  Equations  (29)  and  (31),  modified  as  outlined  in tteae  preceding 

section,  are  therefore  used  to  compute  and  qrev,  resxcootively. 


Before  proceeding  with  the  solution  of  equations  (29)  and  C 3l) , the 

dimensional  character  of  the  structural,  aerodynamic,  integntl.tLng,  and 
interpolating  matrices  involved  in  these  equations  should  betacorj.  side  red. 
The  aerodynamic  and  interpolating  matrices  are  clearly  dimentoC3nJ.es s. 


The  integrating  matrix 


and  the  rolling-moment-arm  matrix  >c 


Lsl  * yjj 


are  proportional  to  the  exposed  semispan  Z.  The  influence-towxc'tion 
matrix  for  the  example  configuration,  which  is  presented  in  tiblole  IV,  is 

written  as  the  product  of  the  quantity  lj—~  and  a nondinwsfcs  ional 
matrix. 


The  aforementioned  independence  of  absolute  dimensions  cm  * no w be 
demonstrated.  Equations  (29)  and  (31)  can  be  written  in  the  fon»xm 


. « 

'‘if*’ 


&r"  ' 4‘P  ‘ :%  /■ 


ft  .414 m.. 
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where  the  multiplication  of  the  various  submatrices  by  the  indicated 
quantities  yields  nondimens ional  results.  Inspection  of  equations  (32) 
and  (33)  shows  that  only  ratios  of  the  dimensions  are  involved,  and  that 
the  quantity 


qcl2 

is  an  important  nondimens ional  parameter  for  this  particular  configuration. 

The  step-by-step  procedure  previously  outlined  has  been  followed 
for  the  example  configuration  and  the  results  are  included  in  figure  3* 
Since  the  dynamic  pressure  in  this  case  is  essentially  an  altitude  param- 
eter, the  results  in  figure  3 are  shown  plotted  against  the  pressure  ratio 
Ph/P0,  where  Pq  is  the  standard  pressure  at  sea  level,  as  obtained  from 

reference  5>  and  Pjj  is  the  free-stream  static  pressure  at  altitude  and 
is  related  to  q by  the  equation 


in  which  7 is  the  ratio  of  specific  heats  of  air 
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Simplified  Methods 

Two  simplified  methods  of  analysis  have  also  been  used  to  solve 
for  qrev  and  0 for  the  example.  They  differ  from  the  method  pre- 
viously described  only  in  that  the  airloads  are  determined  by  simplified 
means.  In  the  first  method,  lifting-surface  theory  is  modified  in  a 
manner  similar  to  that  used  in  reference  6,  wherein  it  is  assumed  that 
the  chordwise  center  of  pressure  due  to  structural  deformation  and  roll 
coincides  with  the  elastic  axis  of  the  wing.  Thus,  the  total  twisting 
moment  about  the  elastic  axis  is  equal  to  the  twisting  moment  due  to 
only  the  aileron,  and  equation  (15)  becomes 


de  = qcf  & rz 
dy  3 J 0 dy 


Pc^fa)  dTl 


(54) 


The  rolling-moment  equation,  equation  (12),  is  unchanged.  A solution 
to  the  problem  has  been  obtained  by  using  matrices  in  a manner  similar 
to  that  employed  previously.  In  this  case,  however,  de/dy  is  given 

explicitly  by  equation  (34)  and,  therefore,  can  calculated 

directly,  the  solution  of  a set  of  simultaneous  equations  being 
unnecessary. 

In  the  second  method  of  analysis,  the  aerodynamic  terms  are  derived 
on  the  basis  of  two-dimensional  (strip)  theory.  For  the  case  wherein 
the  elastic  axis  lies  on  the  midchord,  the  center  of  pressure  due  to 
structural  deformation  and  roll  coincides  with  the  elastic  axis  of  the 
wing  and  the  only  twisting  moment  about  the  elastic  axis  is  that  produced 
by  the  aileron.  Therefore,  equation  (34)  expresses  d0/dy  exactly,  and 
the  rolling  rate  obtained  by  solving  equations  (12)  and  (34)  Is  exact. 
Because  of  the  simplicity  of  the  strip  theory,  these  calculations  can  be 
performed  analytically. 

The  results  obtained  by  using  these  methods,  termed,  respectively, 
"modified  lifting-surface  theory"  and  "strip  theory,"  are  also  shown  in 
figure  3 for  comparison  with  the  results  obtained  by  using  lifting-surface 
theory. 


RESULTS  AND  DISCUSSION 


The  results  shown  in  figure  3 indicate  that  the  variation  of  rolling 
effectiveness  0 with  the  pressure  ratio  P^/Po  for  Practical  pur- 
poses, linear.  This  linearity  suggests  the  possibility  that,  for  configu- 
ations  of  the  type  considered,  only  the  computation  of  pressure  ratio  at 


-s* 


s4i-v*v  S-T  -4*  ' t 
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reversal  (Ph/po)r  v would  required;  the  rolling  effectiveness  for 

smaller  values  of  P^/Pq  could  then  be  obtained  by  linear  interpolation. 

In  this  way  the  computation  time  could  be  greatly  reduced,  the  calcu- 
lation of  (^h/^o)rev  ^e^n8  a relatively  simple  process.  Some  care 

should  be  taken,  however,  in  making  use  of  this  linearity,  since  it 
appears  to  depend  on  the  proximity  of  the  elastic  axis  to  the  chordwise 
center  of  pressure  due  to  angle  of  attack.  For  wings  wherein  the  elastic 
axis  is  distant  from  this  center  of  pressure,  considerable  curvature  of 
the  rolling-effectiveness  curves  could  result.  At  the  other  extreme, 
the  results  in  figure  3 obtained  by  using  the  two  simplified  methods  are 
exactly  linear,  since,  in  both  cases,  the  elastic  axis  and  the  center  of 
pressure  are  coincident. 

An  additional  consequence  of  this  linearity  is  that  the  accuracy 
with  which  a particular  method  predicts  (^h/^o)rev  *s  a ^rec^  measure 

of  its  ability  to  predict  rolling  effectiveness.  A comparison  of  the 
values  of  (^h/^o)rev  °^^a^ne<^  by  the  two  simplified  methods  with  those 

obtained  by  the  method  presented  in  this  paper  is  therefore  shown  in  fig- 
ure 4,  wherein  (Ph/po)rev  is  against  Mach  number.  Although  the 

results  show  very  little  difference  in  the  values  of  v as 

obtained  by  the  three  methods  at  high  values  of  M,  considerable  differ- 
ence exists  at  low  values  of  M.  The  results  obtained  by  the  use  of  the 
modified  lifting-surface  theory  are  consistently  unconservative;  the 
results  obtained  by  the  use  of  strip  theory  are  consistently  conserva- 
tive. Modified  lifting-surface  theory  neglects  the  twisting  moments 
arising  from  the  twist  of  the  wing;  only  the  torques  caused  by  aileron 
deflection  are  considered.  Neglect  of  the  twisting  moments  due  to  angle- 
of -attack  changes  evidently  reduces  the  resultant  angle  of  twist  and 
therefore  reduces  the  adverse  rolling  moments  caused  by  the  deformation. 
When  strip  theory  is  used,  the  elastic  twist  is  again  lower  than  that 
obtained  by  lifting-surface  theory  because,  again,  only  the  twisting 
moments  caused  by  the  aileron  deflection  are  present.  Here,  however, 
the  absence  of  the  finite -span  effects  actually  results  in  a greater 
adverse  rolling  moment,  even  though  the  twist  causing  the  adverse  rolling 
moment  is  in  itself  smaller. 


Further  mention  should  be  made  of  the  behavior  at  large  Mach  numbers. 
Actually  all  three  methods  should  yield  the  same  results  as  M approaches 
infinity  because,  as  M increases  indefinitely,  lifting-surface  theory 
approaches  strip  theory.  As  an  illustration  of  this  fact,  the  rate  of 


roll  for  a rigid  wing 


is  shown  as  a function  of  M in  figure  5. 


obtained  by  using  lifting-surface  theory  rapidly 


The  values  of 
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approach  the  values  obtained  by  using  strip  theory  as  the  Mach  number 
Increases.  Thus,  at  high  values  of  M the  use  of  strip  theory  should 
yield  aeroelastic  solutions  of  an  accuracy  comparable  to  that  obtained 
by  using  lifting-surface  theory. 

In  actual  application  it  may  be  convenient  to  have  rolling  effec- 
tiveness given  as  a plot  of  pressure  ratio  against  Mach  number  for  con- 
stant values  of  rolling  effectiveness.  A graph  of  this  type  may  be  made 
by  cross -plotting  the  information  contained  in  figure  3 and  is  included 
in  figure  6.  Also  included  in  figure  6 is  an  additional  ordinate  that 
gives  the  altitude  as  obtained  from  the  standard-atmosphere  table  in 
reference  5. 

The  rolling  effectiveness  0 at  any  time  during  a particular  flight 
may  be  determined  as  a function  of  Mach  number  if  a history  of  the  flight 
is  known  in  the  form  of  a plot  of  standard  altitude  against  Mach  number. 
For  example,  consider  two  constant -altitude  flights,  one  at  30,000  feet 
and  the  other  at  20,000  feet.  The  resulting  variation  of  the  rolling 
effectiveness  of  the  example  configuration  with  Mach  number  is  shown  in 
figure  7« 

The  indicial-solution  (aerodynamic-influence-function)  approach  used 
herein  in  the  calculation  of  the  aerodynamic  loads  exhibits  considerable 
promise  of  being  applicable  to  plan  forms  other  than  rectangular.  The 
delta  wing  and  the  low-aspect-ratio  swept  wing  at  supersonic  speeds  can 
be  handled  in  the  same  way  as  the  rectangular  wing;  care  should  be  taken 
in  these  cases,  however,  to  account  for  chordwise  bending  where  necessary. 
Other  static  aeroelastic  problems  such  as  center-of -pressure  shift  (which 
has  been  considered  by  Frick  and  Chubb  in  ref.  1 for  high-aspect-ratio 
swept  wings)  and  load  distribution  seem  to  be  amenable  to  analysis  by  the 
methods  contained  herein.  It  might  even  be  possible  to  extend  the 
approach  to  take  into  account  unsteady  aerodynamic  effects  and  thereby 
to  obtain  accurate  solutions  to  flutter  problems.  The  calculations  in 
the  last  case  would  undoubtedly  be  arduous  and  the  main  usefulness  of 
the  approach  would  be  to  establish  a basis  for  the  evaluation  of  more 
practical  but  necessarily  less  accurate  solutions  of  the  flutter  problem. 

Problems  involving  configurations  about  which  the  flow  is  not  sub- 
stantially potential  are  generally  not  amenable  to  this  type  of  approach. 
More  specifically,  the  success  of  the  approach  depends  on  the  applica- 
bility of  linearization  to  the  aerodynamics  and  on  the  ability  to  calcu- 
late the  necessary  indicial  load  distributions. 


CONCLUDING  REMARKS 


A method  has  been  developed  for  the  prediction  of  the  aeroelastic 
effects  on  the  roll  of  rectangular  wings  in  supersonic  flow.  The  method 
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is  based  on  the  use  of  influence  functions,  either  analytically  or 
experimentally  obtained,  to  calculate  the  structural  distortion  due  to 
airloads.  The  airloads,  themselves,  are  calculated  on  the  basis  of  line- 
arized lifting-surface  theory  by  superposing  basic  airloads  resulting 
from  elementary  angle -of -at tack  distributions.  The  solution  of  the  aero- 
elastic  equations  has  been  obtained  by  means  of  a numerical  procedure 
suitable  for  use  with  desk-type  calculators. 

Results  for  an  example  configuration  indicate  that  the  variation  of 
rolling  effectiveness  with  the  free-stream  static  pressure  at  a constant 
Mach  number  is  almost  linear;  a good  approximation  may  be  made  by  assuming 
a linear  variation.  Thus,  in  any  other  cases  wherein  this  linearity  can 
be  expected  - that  is,  when  the  elastic  axis  is  near  the  center  of  pres- 
sure due  to  angle  of  attack  - the  calculations  may  be  greatly  simplified 
in  that  it  is  necessary  to  compute  only  the  free-stream  static  pressure 
at  aileron  reversal. 

The  results  obtained  by  using  the  method  of  analysis  presented  in 
this  paper  are  compared  with  the  results  obtained  by  using  simplified 
methods  of  analysis.  Although  aerodynamic  strip  theory  is  valid  at  high 
Mach  numbers,  too  conservative  results  are  obtained  at  low  supersonic 
Mach  numbers  because  of  neglect  of  finite-span  effects.  A modified 
lifting-surface  theory,  in  which  twisting  moments  due  to  structural 
deformation  and  roll  have  been  neglected,  yields  results  which  are 
unconservative  at  law  supersonic  Mach  numbers. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va. , December  2,  1953- 
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APPENDIX  A 


STRUCTURAL  ANALYSIS  OF  THE  EXAMPLE  CONFIGURATION 


The  deformations  of  the  wing  of  the  example  configuration  are  ana- 
lyzed herein  by  using  the  approximate  plate  theory  presented  in  refer- 
ences V and  7«  Ihis  method  expresses  the  deflection  w of  a plate  at 
any  point  (x,y)  in  the  form  w(x,y)  = W(y)  - x 0(y)  where  W(y)  is  the 
deflection  of  the  wing  at  the  midchord  and  8(y)  is  the  twist  of  the 
wing.  (See  fig.  1.)  An  energy  solution  of  the  problem  is  used  wherein 
the  expression  for  the  potential  energy  of  the  plate  is  written  in  terms 
of  the  approximate  deflection  function.  Minimization  of  this  potential- 
energy  expression  yields  two  ordinary  differential  equations  in  W(y) 
and  0(y).  Since  the  deflection  W(y)  has  no  effect  on  the  airload,  it 
is  eliminated  from  the  two  equations*  the  process  yields  a single  equa- 
tion in  e(y).  This  equation  has  been  derived  in  reference  7 (eq.  (A22) 
of  that  paper)  and  for  a rectangular  cantilever  plate  of  constant  thick- 
ness t becomes,  in  the  notation  of  the  present  paper. 


~ - 2(1  - n)De  P-  = - f 1 M(C)  d£ 


12 


dy> 


ay 


(ai) 


with  the  boundary  conditions 


0(0)  = 


d0(O)  d^e(z) 


dy 


= 0 


dy£ 


In  equation  (Al),  D is  the  plate  stiffness: 

Et5 


D = 


12(l  - p2) 


where  E is  Young's  modulus  and  p is  Poisson's  ratio.  As  could  be 
expected  for  this  structure,  the  twists  are  seen  to  be  dependent  solely 
on  the  twisting  moment  about  the  midchord  M(y) . 

The  solution  to  equation  (Al)  for  a general  M(y)  can  be  obtained 
by  superposing  indicial  solutions  found  by  considering  the  moment  to  be 
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concentrated  at  any  spanvise  station  t).  Bius,  let 

M(y)  = 8(y  - n) 

where  5(y  - n)  is  the  Dirac  delta  function. 

Equation  (Al)  now  becomes 

. 2(1  . m)Do  - -i(,  - y)  (A2) 

12  dy3  W 


where  t)  - y)  is  the  unit  step  function: 


i(n  - y)  = 1 

(y  ^ n) 

■"A 

II 

O 

(y  > i) 

8G«(y,Ti) 

The  rate-of-twist  influence  function  r is  the  quantity  required 

oy 

in  order  to  obtain  the  rate  of  twist  by  superposition.  This  quantity  is 

dG^y^) 

_ m2  1 L0h , _ ^ _ r 

\ - 

8y 

Dc5A2  Coah  \ l- 

t 

sinh  *(l  - ^sinh  ?v 

(y  ^ n) 

► (A3) 

dGj^yjTi) 

- vn*  cosh  ^ ' > )(.,  . „ 1) 

(y  > n) 

dy 

Dc5X2  cosh  * ' 11 

where 

■K  = |^24(1  - \i) 
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With  regard  to  this  solution,  it  should  be  noted  that  the  quantity 
121^/dc^A^  can  be  written  in  the  form  where  G is  the  shear 

modulus  of  elasticity  of  the  material. 

Superposition  of  the  rate-of -twist  influence  function  yields 


d9(y)  = rl 

dy  J0 


M(n)  dq 


Values  of  ^G^/dy  for  l/c  = 3/2  and  p.  = 1/3  (which  yields 
^ = 6)  have  been  computed  for  0 ^ j % 1 in  intervals  of  0.2  and 
0 £ j ^ 1 in  intervals  of  0.1;  the  results  are  included  in  matrix  form 

If 

in  table  IV. 
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APPENDIX  B 

DERIVATION  OF  AERODYNAMIC  COEFFICIENTS 


In  this  appendix  the  aerodynamic  loads  necessary  for  computing 
rolling  effectiveness  are  derived.  The  loads  are  found  in  the  form  of 
section  lift  and  moment  coefficients  by  applying  linearized  supersonic 
lifting-surface  theory  to  the  rectangular  wing  with  three  different 
angle-of -attack  distributions:  the  unit-step  angle-of -attack  distri- 

bution shown  in  figure  2(c),  the  angle -of -at tack  distribution  caused  by 
rolling,  and  the  angle -of -at tack  distribution  resulting  from  aileron 
deflection. 

Analyses  of  each  of  these  problems  are  contained  in  the  literature. 
The  unit-step  problem  is  essentially  the  same  as  the  problem  of  finding 
the  loads  due  to  deflection  of  an  outboard  aileron,  which  has  already 
been  solved  (see,  for  example,  ref.  8);  the  rolling  problem  has  been 
treated  by  many  investigators  (see,  for  example,  ref.  9)*  There  does 
not  seem  to  be  any  report,  however,  that  gives  the  desired  coefficients 
in  a form  sufficiently  complete  for  the  purposes  of  this  paper.  For 
this  reason  and  also  for  convenience,  the  necessary  aerodynamic  quan- 
tities are  derived  completely  herein. 

When  linearized  lifting-surface  theory  is  used,  the  lift  per  unit 
area  of  a thin  wing  is  given  by 


4q  dqp(x,y) 

p(x,y)  = — — (Bl) 

V dx 

where  x is  positive  in  the  direction  of  the  airstream.  The  poten- 
tial cp,  evaluated  at  the  surface  of  the  wing,  is 


fOw)  - I ff  : “S  d?  — (K) 

” 3 |/(l  - *)*  - 8*(C  - ,)2 

where  o(£,£)  is  the  local  angle  of  attack  of  the  wing.  In  general,  for 
a rectangular  wing,  the  region  of  integration  S includes  the  entire  area 
on  the  wing  within  the  forward  Mach  cone  from  the  point  (x,y).  However, 

In  order  to  obtain  the  potential  at  a point  near  a wing  tip,  the  proper 
region  of  integration  S is  determined  by  using  Eward's  method  (ref.  10). 


'*?!•'  ■ . 4 
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The  spanvise  loading  and  twisting  moment  about  the  midchord  are 
given  by 

-/2 


pc/z 

L(y)  = / p(x,y)  dx 
J -c/2 

rc/2 

M(y)  = - x p(x,y)  dx 
J-c/2 


(B3) 


If  the  expression  for  p(x,y)  is  substituted  into  equations  (B3),  and 
the  results  are  integrated,  the  lift  and  moment  become: 


L(y)  - ^ <p(|.y) 


M(y)  - -ia 


|»(|,y)-^j%(x,y> 


c/s 


(b4) 


where  the  fact  that  the  potential  is  zero  at  the  leading  edge  has  been 
taken  into  account. 

The  foregoing  equations  can  be  nondimens ionali zed  by  letting 


. Lp 
Be.  = — 

K I qc 

qc^ 

1 . x 

xi  = x + - 
1 2 c 


4 = — 9 
VI 


J »'  . 

\ ■ 

, -V 


V i 


• . ? aAkssL 


*¥  ii 

.if 


• \ 


#4 
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The  resulting  equations  are 


tci(y 1)  - t *(ltyi) 


pcm(yi)  -™ 


\ - f *(xi>yi)  <**i 


(B5) 


where 


*(*,,»!)  - /r 

8 ^|(xi  - lx)2  - m2(yi  - 5)S 


(b6) 


and  m is  the  modified  aspect-ratio  parameter  32/c . 


Derivation  of  pcjg'  and  3cme' 


The  equations  for  the  spanvise  loading  pCjg’  and  the  twisting 
moment  about  the  midchord  Pcmo ' due  to  the  unit-step  angle -of -attack 
distribution  shown  in  figure  2(c)  become,  from  equations  (B5), 


km 


(B7) 


pcme’Cyi^i) 


4m  1 

x 2 


•(^yi»ni) 


4(xi>yi^i)  ^1 


where  the  potential  4>  is  found  for  an  angle -of -at tack  distribution 
defined  by  the  unit  step  function: 


a(yi^i)  =2(n  1 - yi) 


a(yi,ii)  = 1 

a(yi,tji)  = 0 


or 


fri  £ yi) 

0»i  < yi) 
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^3 


Equation  (B 6)  136000168 


dCi  dfti 


/T_ 

E |/(xi  - li)2  - m2(yi  - ^ 


(B8) 


The  region  of  integration  E includes  only  that  portion  of  S which 
lies  within  the  region  of  the  wing  between  the  wing  tip  and  the  position 
of  the  unit  step  in  angle  of  attack  at 


The  plan  form  and  angle -of -attack  distribution  considered  are  those 
shown  in  figure  2(c).  From  this  figure  it  can  be  seen  that  the  position 
of  the  unit  step  is  restricted  to  the  range  0 £ £ 1.  If,  in  addi- 

tion, the  modified  aspect-ratio  parameter  m is  restricted  to  values 
greater  than  ■ ^ — , the  potentials  4 for  the  finite  wing  in  fig- 
ure 2(c)  are  exactly  the  same  as  the  potentials  on  the  finite  wing  por- 
tion of  the  following  semi -infinite  plan  form  which  has  been  obtained  by 
allowing  the  left  tip  and  the  trailing  edge  of  the  finite  wing  to 
approach  infinity: 


YI 

'n  1^- 

nnof  ir 

rz(o) 

s' 

^^rz(b)  ^ 

* ^ ^ 

'''  3C(o) 

Thus,  attention  is  confined  to  the  derivation  of  the  potentials  on  this 
semi-infinite  plan  form.  With  the  potentials  known,  the  loading  coeffi- 
cients on  the  finite  wing  are  given  by  equations  (B7). 

The  Mach  lines  shown  in  the  above  sketch  bound  a number  of  distinct 
regions.  These  regions  are  significant  because  the  area  of  integration 
in  equation  (b8)  takes  a different  geometric  form  in  each;  it  is  to  be 
expected,  therefore,  that  the  potential  4 is  given  by  a different 
equation  in  each  region.  The  proper  area  of  integration  L (shown 
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In  region  III(A) : 


>111 (A) 


tan" 


In  region  III(b)  : 


IX1  + ml 

^1  - yi) 

IX1  - “1 

(’ll  - yi) 

.-1  X1 


|xi  + m(Tii  - yi) 


X 


1 1 


•m(B)  = 5 


+ m 

(’ll  - yi) 

x!  - m 

(’ll  - yi) 

• ! 
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Notice  that  the  only  difference  in  the  equations  for  *m(A)  and 

Is  that  the  radical  in  the  inverse -hyperbolic -tangent  term  in 
the  equation  for  #III(A)  is  the  reciprocal  of  the  radical  in  the 
corresponding  term  in  the  equation  for  ®ui(b)*  Thus,  the  following 
form  holds  for  both  regions: 


2 1 X1  + “Hi  - yi)  1 

*111  = 5 X1  tan  Jxx  - yi)  + “(^l  “ yi)tanh  kl 


where 


. xi  - nOi  - yi) 

Pi  + mp!  - yi) 

- Jxi + m(ni  ~ ylT 

' !X1  - “(Hi  - yi) 


(yi  S \) 


(*i  * \) 


A similar  situation  occurs  in  the  derivation  of  ®jv(a)  ancl  ®IV(B) 
and  of  and  ^V(B)  * Therefore,  the  equations  for  G’jy  and 

may  be  derived  for  either  the  subregion  (A)  or  (B)  and  applied  through- 
out the  entire  region. 

In  region  IV: 


y.)' 


n 

2 1 — i -ilxi + “di  - yi)  .ilxi  • *wi  .1. 

dlKiK  - ^1  + X1  - m(ni~-  yi)  ' *“  + a(ni  ' yi)tanh  ^ 

w - 


48 


NACA  TN  3067 
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Since  the  spanwlse  loading  3cig'  Is  determined  by  the  value  of  $ 
at  the  trailing  edge  (see  eq.  (B7)),  the  equation  for  3c jg'  is 

dependent  only  on  the  region  that  the  trailing  edge  is  in.  For  example, 
if  for  a given  value  of  yj_  the  trailing  edge  is  in  region  III,  the 
value  of  is 

|fcl0'(yl>’ll2]iii  = ~ ®IllOrl»1»lli) 


and  the  particular  form  to  use  for  3c is  determined  by  the  region 
in  which  the  trailing  edge  falls. 
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^9 


Ohe  determination  of  requires  evaluation  of  the  integral 

of  the  potential  along  the  chord  of  the  wing.  However,  the 

hanS^o  bf  the+fXpreB8i0ns  for  the  Potentials  valid  within  each  region 
happen  to  be  continuous  across  the  boundaries  of  each  region  Furthfv 

^flnlte  ‘”*•"1  of  the  potential  th“  le 
neighborhood  of  the  leading  edge  Is  sero  at  the  leadlna  edae  eC  th„ 
Integral  of  the  potential  along  the  chord  may  be  found  by  evaluating  at 
the  trailing  edge  the  Indefinite  Integral  of  th'  potential  ^^.1^ 

trailing  edge  falls  successively  in  each  of  the  regions  I to  VI  the 
resulting  equations  for  go.  are,  therefore,  ’ 


K>V 0 


’)n  * T.{1  - ^i)5/2^I 


')m  ■ ; "("1  - »i) 


tanh 


-K  V1  - m - yi)! 

K1  " X 


(SV)rv  - H 


m(i,i  - yi) 


tanh"1^  - £ : m2(^  ' yl £.  + l (B10) 


f(l  - 


Jr 


J 


')  ■- 
)xr  Jt 


(X  " "T^)ft^i  + 0(11!  - y1)tanh“1kg 


’)  - 0 

'VI 
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i 


1 - mfrj!  - yx) 


1 + m(n!  - yx) 


J& 

hi 


jl  + m 


(’ll  - yi) 


(^1  - yi) 

As  before,  the  particular  form  to  use  for 
region  in  which  the  trailing  edge  falls. 


(*1  * 1l) 

(*i  i ii) 

' is  determined  by  the 


Equations  (B9)  and  (BIO)  may  now  be  used  to  obtain  the  spanwise 
loading  and  twisting  moment  about  the  midchord  due  to  a unit-step  angle- 
of -attack  distribution  for  particular  values  of  and  m.  In  summa- 
rizing the  information,  it  is  convenient  to  separate  m and  q^  into 
various  ranges.  Sketches  of  the  Mach  lines  and  summaries  of  the  par- 
ticular forms  to  use  in  each  range  follow: 

For  m ^ 2 and  0 ^ q^  ^ i; 


for  region  III  when  j S y1  = ~ + 


*1 


m 


^1 


m + ^1  - yl  - 2(1  + a) 


for  region  VI  when 
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Use  the  coefficients  < 


for  region  II  when  0 ^ y,  £ — 

x m 

for  region  I when  ^ ^ y,  * n,  - i 

mi.  x m 

for  region  III  when  ri-  - i < y < n,  + i 

x m x *x  m 

for  region  VI  when  ^1  + 5 - yi  - 2(1  + a) 


For 


— - — % m ^ 2 and  0 ^ Hi  - — : 
1 + 2a  1 " o 


Use  the  coefficients  < 


for  region  V when  0 y,  < i - ti. 

1 in  1 

for  region  IV  when  — - n,  S y,  $ ^ 

in  x x in 

for  region  III  when  - ^ y,  ^ — + rj, 

in  x in  x 

for  region  VI  when  i + tj1  % % 2(l  + a) 
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For  — - — S m S 2 and  - £ tj  < Is 
1 + 2a  ml 


Use  the  coefficients  •< 


for  region  II  when  0 £ - i 

for  region  IV  when  ^ - yl  - m 

for  region  III  when  ^ ^ yl  “ ^1  + m 

for  region  VI  when  i % y1  ^ 2(1  + 


Derivation  of  pc^  and  (Jc^ 

The  local  angle  of  attack  due  to  a unit  rate  of  roll  (ph/2V 
is  given  by 


- -(l  - ^77) 


The  potential  is  therefore  given  by 


®(xl^l)  = -//  I1 
3 


?1 ) 

dCi  dlj. 

1 + ■‘Act- 

. *.\2  -2/v  _ r. \2 

a) 

• 1) 
(Bll) 

(B12) 


f 
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In  this  section,  the  potential  $ is  derived  for  a finite  rectan- 
gular wing  with  a nondimensional  span  of  2(1  + a) . In  order  to  simplify 
the  derivation,  it  is  assumed  that  no  point  on  the  wing  is  influenced  hy 

hoth  wing  tips;  this  assumption  is  valid  provided  that  m £ —-i — . 

1 + a 

TOe  Mach  lines  originating  at  the  intersections  of  the  leading  edge 
and  the  tips  divide  the  wing  into  two  distinct  regions  (see  the  following 
sketches)  and  the  potential  ♦ for  these  two  regions  is  given  toy  equa- 
tion (B12)  as  follows: 


In  region  I: 


In  region  II: 
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®II 


2 1 
m 1 + a 


- T ' I t)^iK  • "1,1  + “it1 


°yi 

xi  “ “Vl 


2 a 
m 1 + a 


X1  - "Vi 


+ x. 


°yi  \ 

xi  - "Vl  / 


The  equations  for  the  s panwise  loading  3c ^ and  the  twisting  moment 

about  the  midchord  3cmp  are  found  in  a manner  similar  to  that  used  to 

determine  the  equations  for  3c  7 » and  3cm  lkus, 

*0  0 


The  equations  for  3^,  obtained  in  a manner  similar  to  that  used  pre- 
viously, are 


(%)x  * 0 


8 1 
x 1 + a 


y W) 


8 a 
x 1 + a 
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In  order  to  conform  to  the  assumption  made  in  the  aerodynamic 
analysis,  the  airloads  over  the  rigid  plate  in  the  center  section  of 
the  wing  are  to  be  neglected.  Therefore,  the  limit  of  the  modified 

aspect-ratio  parameter  can  be  lowered  from  m « - — to  m £ 

1 + a l+2a 

Summaries  of  the  particular  forms  to  use  and  illustrative  sketches 
follow: 

For  m £ 1: 


Use  the  coefficients  < 


for  region  II  when  0 % y^  % i 
for  region  I when  g % y^  $ 1 


Use  the  coefficients  for  region  II  when  0 £ y1  £ 1 


NACA  TN  3067 


It  should  he  noted  that  each  expression  for  0c  and  pc^  can 

he  separated  into  two  parts  - one  with  the  coefficient  — i — and  the 

1 + a 

other  with  the  coefficient  — each  of  which  is  independent  of  a. 

1 + a 

The  first  part  is  the  result  that  would  be  obtained  on  the  right  wing  if 
the  wing  were  rolling  about  the  x-axis  (see  fig.  2(a));  it  is  therefore 
designated  0c  ^ or  pc^  . The  second  part  can  be  shown  to  be  the 

section  lift  or  moment  coefficient  which  results  from  a uniform  angle 


therefore  designated 

^la  or  Pc% 

00^  can  be  written 

as 

a Be 

C*p  1 + a PC*p0 

1 + a *a 

. a Bc 

PC“P  " 1 + a PCmP0 

1 + a 

where  the  coefficients  in  the  various  regions  are: 

(‘"'Jpjj  ’ J*(1  - n) 
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Derivation  of  3c  and  30^ 


In  this  section,  the  aerodynamic  loads  due  to  aileron  deflection 
are  derived  for  a trailing -edge  aileron  of  constant  chord  ca  and 
span  ba.  The  spanvise  location  of  the  inboard  end  of  the  aileron  is 
left  general;  the  outboard  end  of  the  aileron  coincides  with  the  wing 
tip.  It  is  assumed  that  there  are  no  gaps  between  the  aileron  and  the 
wing.  Furthermore,  the  portion  of  the  wing  forward  of  the  hinge  line 
may  be  assumed  to  be  absent  because  the  wing  is  considered  to  be  at  zero 
angle  of  attack,  the  aileron  hinge  line  is  unswept,  and  only  supersonic 
speeds  are  considered.  The  resulting  configuration  is  exactly  the  same 
as  that  encountered  in  the  derivation  of  3c jq  and  0cmg.  In  addition, 

since,  for  most  reasonable  aileron  configurations,  the  pressures  on  the 
right  wing  are  unaffected  by  the  deflection  of  the  left  aileron,  the 
total  spanwise  loading  3c  jg  and  the  total  twisting  moment  about  the 

midchord  3cm^  can  be  found  by  properly  substituting  the  aileron  dimen- 
sions for  the  wing  dimensions  in  the  results  already  obtained  for  0c jg' 
and  0cmg'.  In  order  that  the  foregoing  be  true,  the  restriction 
^ C&1 

m s — = is  necessary. 

1 + 2a  - b&1 


If  ca  is  substituted  for  c and  b^  = ba/z  for  in  the 

expressions  for  the  loads  due  to  structural  deformation,  the  load  due 
to  a unit  aileron  deflection  will  result.  Several  more  steps  are  neces- 
sary, however,  in  order  to  produce  the  desired  results  3c^B  and  3cmg' 

In  the  first  place,  the  reduced  aspect-ratio  parameter  which  was  3Z/c 
becomes  0Z/ca  = m/cai,  where  cai  = ca/c.  In  addition  to  the  substi- 
tutions mentioned  in  the  preceding  paragraph,  therefore,  the  substitu- 
tion of  nyca^  for  m must  be  made,  both  in  the  expressions  for  0c jg* 

and  00*  and  in  the  expressions  which  define  the  limits  of  applies- 
0 


bility  of  these  terms.  Furthermore,  in  the  nondimensionalization  of  the 
moment  expression  the  quantity  c is  involved.  In  order  to  preserve  the 


nondlmens ionali z ing  coefficient  qc/0  for  the  aileron  loads  and  liC2/0 
for  the  aileron  moments,  the  egressions  for  0cZq  ' (yi^ay)  and 

("Vfs'i-v)  must  be  multiplied  by  cai  and  c&12,  respectively. 


Lastly,  since  the  moments  are  taken  about  the  midchord  in  the  case  of 
pressures  due  to  the  structural  deformation,  direct  substitution  yields 
the  moment  about  the  aileron  midchord;  thus,  the  moment  must  be  trans- 
ferred back  to  the  wing  midchord. 


cai 

for  region  II  when  0 % 


Use  the  coefficients  ■% 


for  region  I when 


m 


*1  m 


^a-i  ^ ^ ca-i 

-=±S  yx  = ba,  - -=* 
al  < v < 


for  region  III  when  ba 3 y,  = 1 

1 m 


The  expressions  for  the  coefficients  in  the  various  ranges  are: 
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TABLE  I.  - PARTIAL  SECTION  LIFT  AND  M3MWT  COEFFICIENTS  Bcle'  AND  Sc^1 
FOR  A IHIT-STEP  ANCLE  OF  ATTACK 


(a)  $j/e  . 5/7 


y/i 

For  tj/l 

equal  to  - 

_ 0 

0.1 

0.2 

0.3 

o.lf 



0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

0cia’(y/i»V!) 

-1.3 

0.025339 

0 

0 

0 

0 

0 

0 

0 

0 

— 

0 

0 

-1.2 

.067298 

.023359 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.1 

.126173 

.067298 

.023339 

0 

0 

0 

0 

0 

0 

0 

0 

-1.0 

.198493 

.126173 

.067298 

.023539 

0 

0 

0 

0 

0 

0 

0 

-.9 

.283854 

.198495 

.126173 

.067298 

.023339 

0 

0 

0 

0 

0 

0 

-.8 

.382447 

.283854 

.198493 

.126173 

.067298 

.023339 

0 

0 

0 

0 

0 

-7 

.494931 

.382447 

.283854 

.198493 

.126173 

.067298 

.023339 

0 

0 

0 

0 

-.6 

.622469 

.1,94951 

.382447 

.283854 

.198493 

.126173 

.067298 

.023339 

0 

0 

0 

-.5 

.7^5859 

.622469 

.494931 

.382447 

.283854 

.198493 

.126173 

.067298 

.023359 

0 

0 

-.1* 

.930849 

.766859 

.622469 

.494931 

.382447 

.283854 

.198493 

.126173 

.067298 

.023339 

0 

-3 

1.04930 

.897717 

.733726 

.589337 

.461799 

.349315 

.250723 

. 165362 

.093041 

.034167 

0 

-.2 

1.24247 

1.022886 

.834923 

.670937 

.526543 

.399010 

.286521 

.187932 

.102572 

.035611 

0 

-.1 

1.42537 

1.157726 

.938137 

.750178 

.586189 

.441798 

.314260 

.201778 

. 107498 

.037253 

0 

0 

1.71432 

1.320558 

1.052711 

.833126 

.645163 

.481175 

.336785 

.213062 

.113216 

.039149 

0 

.1 

1.98515 

I.589087 

1.195130 

.927484 

.707896 

.519936 

.359491 

.226504 

.119961 

.041365 

0 

.2 

2.10430 

1.836762 

1.442695 

1.048740 

.781092 

.564900 

. 387676 

.242903 

.128081 

.044005 

0 

.3 

2.15443 

1.934844 

1.667305 

1.273240 

.882613 

.625377 

.424091 

.263557 

.138125 

.047229 

0 

.u 

2.14681 

1.958854 

1.739267 

1.474968 

1.091348 

.715331 

.474109 

.290757 

.151004 

.051285 

0 

.5 

2.08438 

1.920391 

1.735864 

1.526792 

1.277190 

.909457 

.550858 

.329142 

.168392 

.056609 

0 

.6 

1.96398 

1.823206 

1.670192 

1.500869 

1.308059 

1.076031 

.727565 

.390548 

.193838 

.064041 

0 

.7 

1.77946 

1.663767 

1.539261 

1.403506 

1.252644 

1.079785 

.869634 

.545674 

.237054 

.075502 

0 

.8 

1.51367 

1.423440 

1.327050 

1.223027 

1.109179 

.981991 

.835096 

.654029 

.363783 

.096919 

0 

.9 

1.11123 

1.049983 

.904929 

.915250 

.839798 

.756835 

.663525 

.554592 

.417548 

*-* 

CD 

« 

M 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Sc^'Cy/l.n/l) 

-1.3 

-0.010999 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.2 

-.029761 

-.010999 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1.1 

-.052085 

-.029761 

-.010999 

0 

0 

0 

0 

0 

0 

0 

0 

-1.0 

-.076012 

-.052085 

-.029764 

-.010999 

0 

0 

0 

0 

0 

0 

0 

-.9 

-IOOO87 

-.076012 

-.052005 

-.029761 

-.010999 

0 

0 

0 

0 

0 

0 

-.8 

-.123018 

-.100087 

-.076012 

-.052085 

-.029761 

-.010999 

0 

0 

0 

0 

0 

-.7 

-.143537 

-.125018 

-.100087 

-.076012 

-.052085 

-.029761 

-.010999 

0 

0 

0 

0 

-.6 

-.160289 

-.143537 

-.123018 

-.100087 

-.076012 

-.052085 

-. 029761 

-.010999 

0 

0 

0 

-.5 

-.171703 

-.160289 

-.143537 

-.123018 

-.100087 

-.076012 

-.05208? 

-.  029761 

-.010999 

0 

0 

-.It 

-.175804 

-.171703 

-.160289 

-.143537 

-.123018 

-,100087 

-.076012 

-. 052085 

-.029761 

-.010999 

0 

-.3 

-.154246 

-.160190 

-.156087 

-.144674 

-.127922 

-. 107401 

-.084472 

-.060397 

-.036470 

-.014145 

0 

-.2 

-.107098 

-.127428 

-.133372 

-.129270 

-.117857 

-.101104 

-,080584 

-.057654 

-.033579 

-.012198 

0 

-.1 

-.031475 

-.074895 

-.095225 

-. 101170 

-.097067 

-.085653 

-.068902 

-.048381 

-.027503 

-.010100 

0 

0 

.109560 

.003450 

-.039971 

-.060299 

-.066244 

-.062141 

-.050727 

-.  035792 

-.020862 

-.007815 

0 

.1 

.251367 

.145258 

.039147 

-.004272 

-.024602 

-.030546 

-.028130 

-.021767 

-.013514 

-.005301 

0 

.2 

.329556 

.286136 

.180025 

.073915 

.050495 

.008549 

.002149 

.005845 

.005251 

.002495 

0 

.3 

.382066 

.361737 

.318317 

.212207 

.104509 

.055480 

.028715 

.012702 

.004245 

.000698 

0 

a 

.415838 

.409894 

.389565 

.344554 

.233860 

.120461 

.067324 

.035155 

.015493 

.004426 

0 

.5 

.433166 

.437268 

.429691 

.404701 

.353949 

.238191 

.120487 

.063932 

.029415 

.008936 

0 

.6 

.434189 

.443882 

.443121 

.429608 

.399434 

.344338 

.225199 

.105261 

.047922 

.014697 

0 

• 7 

.415320 

.426825 

.430171 

.425910 

.405836 

.372163 

.314800 

.194884 

.076094 

.022748 

0 

.8 

.369481 

.380991 

.386357 

.384655 

.374563 

.354045 

.319554 

.265366 

.147246 

f 

w 

0 

0 

• 9 

.281716 

.291131 

.296477 

.297212 

.292601 

.281579 

.262471 

.232288 

. 184276 

.082284 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6k 
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TABLE  I.-  PARTIAL  SECTION  LOT  AND  M3KENT  COEFTICIENTS  Scj9'  AMD  00^' 


FOR  A UHT-STH*  ANGLE  Or  ATTACK  - Continued 


(b)  01/e  - 1 


y/t 

For  tj/z 

equal  to  - 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

-0.9 

0.038956 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.8 

.115297 

.058956 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.7 

.214524 

.113297 

.038956 

0 

0 

0 

0 

0 

0 

0 

0 

-.6 

.341391 

.214524 

.113297 

.038956 

0 

0 

0 

0 

0 

0 

0 

-.5 

.1*94934 

.341391 

.214524 

.113297 

.058956 

0 

0 

0 

0 

0 

0 

-.4 

.678079 

.494934 

.341591 

.214524 

■ 113297 

.038956 

0 

0 

0 

0 

0 

-.3 

.896507 

.678079 

.494934 

.341591 

.214524 

. 113297 

.038956 

0 

0 

0 

0 

-.2 

1.159058 

.896307 

.678079 

.494934 

.541391 

.214524 

.113297 

.038956 

0 

0 

0 

-.1 

1.491553 

1.159^8 

.896307 

.678079 

.494934 

•341391 

.214524 

.115297 

.038956 

0 

0 

0 

2.000000 

1.491351 

1.159860 

.896307 

.678076 

.494934 

.341389 

.214524 

. 113298 

.038956 

0 

.1 

2.453263 

1.944614 

1.435965 

1.104472 

.840921 

.622691 

.439545 

.286005 

.159137 

.057904 

0 

.2 

2.678065 

2.546570 

1.837924 

1.329272 

.997784 

.734229 

.516001 

.532855 

.179313 

.061607 

0 

.3 

2.794595 

2.531045 

2.199552 

1.690903 

1.182254 

.850761 

.587210 

.368980 

.193375 

.066116 

0 

.4 

2.825815 

2.607584 

2.544034 

2.012541 

1. 503895 

.995243 

.663753 

.407062 

.211406 

.071799 

0 

.5 

2.778511 

2.595163 

2.376932 

2.113382 

1.781889 

1.273240 

.771201 

.460798 

.235749 

.079250 

0 

.6 

2.649749 

2.496209 

2.313064 

2.094830 

1.831283 

1.506444 

1.018592 

.546767 

.271373 

.089656 

0 

.7 

2.428465 

2.301597 

2.148057 

1.964909 

1.753704 

1.511699 

1.217489 

.763944 

.331875 

.105701 

0 

.8 

2.085961 

1.984756 

1.857870 

1.712237 

1.552851 

1.374788 

1.169134 

.915640 

.509296 

.135688 

0 

.9 

1.544520 

1.469978 

1.378901 

1.281350 

1.175722 

1.059575 

.928935 

.776427 

.584570 

.254648 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^'(y/j.n/J) 

-0.9 

-0.017908 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.8 

-.047458 

-.017808 

0 

0 

0 

0 

0 

0 

0 

0 

0 

— 7 

-.080858 

-.047458 

-.017908 

0 

0 

0 

0 

0 

0 

0 

0 

-.6 

-.114061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

0 

0 

0 

-5 

-.145536 

-,u4o6i 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

0 

0 

-.4 

-.165594 

-.143536 

-.114061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

0 

-.3 

-.175684 

-.165594 

-143536 

-.  114 061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

-.2 

-.167131 

-.175684 

-I65594 

-.143536 

-.114061 

-.080858 

-.047438 

-.017908 

0 

0 

0 

-.1 

-.127212 

-.167131 

-.175684 

-.165594 

-.145536 

-.  H4o6l 

-.080858 

-.047438 

-.017908 

0 

0 

0 

0 

-.127212 

-.167131 

-.175684 

-.165594 

-.143556 

-.114062 

-.080858 

-.047438 

-.017908 

0 

.1 

-.152677 

.025464 

-.101748 

-.141666 

-.150219 

-.140129 

-. 118071 

-.088596 

-.055393 

-.021977 

0 

.2 

•235037 

.195119 

.067905 

-.059306 

-.099224 

-.107778 

-.097688 

-.075629 

-.046156 

-.017213 

0 

• 3 

.292379 

.283825 

.243907 

.116694 

-.010518 

-.050436 

-.058989 

-.048899 

-.030353 

-.011857 

0 

.4 

.331928 

.542019 

•333467 

.293547 

.166356 

.039122 

-.000794 

-.012545 

-.011913 

-.005684 

0 

.5 

.3557*3 

.377801 

.387891 

•379338 

.339419 

.212207 

.081914 

.033149 

.010505 

.001664 

0 

.6 

.363565 

•393039 

.415098 

.4251&7 

.416635 

.373614 

.237672 

.096958 

.039652 

.010875 

0 

• 7 

.353146 

.386348 

.415823 

.437881 

.444699 

.427101 

.373510 

.229183 

.082769 

.023446 

0 

.8 

.319064 

.352484 

.385686 

.411475 

.423574 

.418971 

.392499 

.333068 

.186742 

.045835 

0 

• 9 

.247091 

.276623 

.305315 

.326265 

.338357 

.339983 

.328655 

.300000 

.244267 

.110347 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

TABLE  I,  - PARTIAL  SECTION  LIFT  AND  M5MENT  COEFFICIENTS  3c ,e'  AND  Bc^ 


FOR  A UNIT-STEP  AN3LE  OF  A1TACK  - Continued 


(e)  B!/c  - 4/3 


y/i 

For  f| Jl 

equal  to  - 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

9ci8'(y/i»VO 

-0.7 

0.021016 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.6 

.113290 

.021016 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.5 

.253952 

.113298 

.021016 

0 

0 

0 

0 

0 

0 

0 

0 

-.4 

.440646 

.253952 

.113298 

.021016 

0 

0 

0 

0 

0 

0 

0 

-.5 

.678076 

.440646 

.255952 

.113298 

1 

8 

0 

0 

0 

0 

0 

0 

-.2 

.978405 

.678076 

.440646 

.253952 

.113298 

.021016 

0 

0 

0 

0 

0 

-.1 

1.370753 

.978405 

.678076 

.440646 

.253952 

.113298 

. 021016 

0 

0 

0 

0 

0 

2.000000 

1.570755 

.978405 

.678076 

.440646 

.255952 

.113298 

.021016 

0 

0 

0 

.1 

2.629246 

2.000000 

1.370753 

.978405 

.678076 

.440646 

.253952 

.113298 

.021016 

0 

0 

.2 

3.021605 

2.629246 

2.000000 

1.370753 

.978395 

.678076 

.440646 

.253952 

.113298 

.021016 

0 

.3 

3.292095 

2.991777 

2.599419 

1.970174 

1.340925 

.948567 

.648249 

.410818 

.224124 

.083470 

0 

.4 

3.397278 

3.159848 

2.859530 

2.467172 

1.837926 

1.208678 

. 816320 

.516002 

.278571 

.095730 

0 

.5 

3.379164 

3.192470 

2.955040 

2.654721 

2.262365 

1.633118 

1.003870 

.611522 

.314331 

.IO5669 

0 

.6 

3.242058 

3.101401 

2.914707 

2.677276 

2.376958 

1.984600 

1.355355 

.729023 

.361851 

.119542 

0 

• 7 

2.970118 

2.877841 

2.737187 

2.550492 

2.313062 

2.012737 

1.623319 

1.018592 

.442500 

.140936 

.8 

2.507915 

2.507915 

2.394621 

2.253963 

2.067269 

1.833051 

1.558846 

1.220853 

.679061 

.180917 

0 

.9 

1.817493 

1.817493 

1.817493 

1.704192 

1.567623 

1.412759 

1.238580 

1.035232 

.779423 

.359531 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Bc^’(y/i»n/») 

-0.7 

-0.009946 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.6 

-.047439 

-.009946 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.5 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

0 

0 

0 

0 

-..4 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

0 

0 

0 

-.3 

-.165594 

-.134544 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

0 

0 

-.2 

-.175340 

-.165594 

-.154344 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

0 

-.1 

-.145363 

-.175340 

-.165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

0 

0 

.0 

0 

-.145363 

-.175340 

-.165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

0 

.1 

.145363 

0 

-.145363 

-.175340 

-.165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

0 

.2 

.175340 

.145363 

0 

-.145363 

-.175340 

-.165594 

-.134344 

-.092126 

-.047439 

-.009946 

0 

.179710 

.189456 

.159479 

.141157 

-.131247 

-.161224 

-. 151478 

-.120229 

-.078010 

-.035323 

0 

.4 

.202250 

.233500 

.243246 

.213269 

.067906 

-.077456 

-.107343 

-.097687 

-.066437 

..026062 

0 

• 5 

.225508 

.267723 

.298974 

.308722 

.278743 

.133380 

.011982 

-.041960 

-.035713 

-.014650 

0 

.6 

.245058 

.299745 

.331962 

.363213 

.372959 

.342982 

.197619 

.050866 

.010187 

-.000590 

0 

.7 

.259451 

.296943 

.341631 

.383848 

.415098 

.424844 

.393465 

.237671 

.073396 

.018193 

0 

.8 

.275268 

.275268 

.322707 

.367394 

.409611 

.439329 

.437969 

.388647 

.218809 

.047776 

0 

•9 

.250073 

.250073 

.250073 

.297512 

.340250 

.368951 

.377847 

.360794 

.304348 

.139585 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABI2  Z.- PARTIAL  SECTIOH  LOT  AHD  HOMERT  COEFFICIOnB  0Cj9'  AMD  00^' 


TOR  A URIT-STEP  AWLE  OP  ATTACK  - Continued 


(d)  0l/c  - 2 


y/i 

Por  ij/l  equ»l  to  - 



0 

0.1 

0.2 

0.5 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

Pe,9'(y/i,n/J) 

-o.u 

0.113298 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.3 

.341391 

.115298 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.2 

.678076 

•341391 

.115298 

0 

0 

0 

0 

0 

0 

0 

0 

-.1 

1.159860 

.678076 

.341391 

.115298 

0 

0 

0 

0 

0 

0 

0 

0 

2.000000 

1.159860 

.678076 

.341591 

. 115298 

0 

0 

0 

0 

0 

0 

.1 

2.840139 

2.000000 

1.159860 

.678076 

•341391 

.113298 

0 

0 

0 

0 

0 

.2 

3.321925 

2.840139 

2.000000 

1.159860 

.678076 

.341391 

.113298 

0 

0 

0 

0 

.3 

3.658608 

3.321923 

2.840139 

2.000000 

1.159860 

.678076 

.341391 

.115298 

0 

0 

0 

.1* 

3.886701 

3.658608 

3.321923 

2.840159 

2.000000 

1.159860 

.678076 

.541591 

.113298 

0 

0 

.5 

4.000000 

3.886701 

5.658608 

3.321925 

2.840139 

2.000000 

1.159860 

.678076 

.341391 

.113298 

0 

.6 

3.387921 

3.837921 

5.724625 

5.496550 

3.159844 

2.678061 

1.857921 

.997782 

.515998 

.179313 

0 

• 7 

3.503889 

3.505889 

5.505889 

3.390596 

3.162503 

2.825813 

2.544029 

1.505889 

.665755 

.211409 

0 

.8 

2.99U40 

2.991140 

2.9911^0 

2.991140 

2.877837 

2.649749 

2.313059 

1.851270 

1.018592 

.271373 

0 

• 9 

2.199262 

2.199262 

2.199262 

2.199262 

2.199262 

2.085964 

1.857870 

1.552848 

1.169154 

.509296 

0 

1.0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^'(y/i,n/i) 

-O.U 

-0.047438 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.3 

-.114062 

-.047438 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.2 

-165592 

-.114062 

-.047458 

0 

0 

0 

0 

0 

0 

0 

0 

-.1 

-.167131 

-.165592 

-.114062 

-.047438 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.167131 

-.165595 

-.114062 

-.047458 

0 

0 

0 

0 

0 

0 

.1 

.167131 

0 

-.167151 

-.165592 

-. 114062 

-.047438 

0 

0 

0 

0 

0 

.2 

.165595 

.167131 

0 

-.167151 

-.165565 

-.114062 

-.047438 

0 

0 

0 

0 

.3 

.114062 

.165592 

.167151 

0 

-.167151 

-.165592 

-.114062 

-.047458 

0 

0 

0 

.4 

.047458 

.114062 

.165595 

.167151 

0 

-.  167131 

-.165595 

-.114062 

-.047458 

0 

0 

• 5 

0 

.047438 

.114062 

.165592 

. 167151 

0 

-.167151 

-.165592 

-.114062 

-.047438 

0 

.6 

.067907 

.067907 

.115545 

.181969 

.233502 

.235055 

.067907 

-O99223 

-.097688 

-.046154 

0 

.7 

.166556 

.166536 

.166556 

.215774 

.280398 

.331928 

.333467 

. 166336 

-.000795 

-.011912 

0 

.8 

.249504 

.249504 

.249504 

.249504 

.296942 

.363566 

.415099 

.416632 

.237672 

.059653 

0 

.9 

.271625 

.271625 

.271625 

.271625 

.271625 

.319064 

.385685 

.423574 

•392499 

.186741 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABLE  I.  • PARTIAL  SECTION  LOT  AMD  MOMENT  COEFFICIENTS  Bc^1  AND  So^' 
FOR  A UNIT-STEP  AN3LE  OF  ATTACK  • Concluded 

(e)  Pl/c  - 4 


y/» 

For  tj/l  equal  to  - 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

-0.2 

O 

B 

CD 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.1 

.678075 

.113298 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2.000000 

.678075 

.113298 

0 

0 

0 

0 

0 

0 

0 

0 

.1 

3.321925 

2.000000 

.678075 

.113298 

0 

0 

0 

0 

0 

0 

0 

.2 

5.886711 

3.521923 

2.000000 

.678075 

.113298 

0 

0 

0 

0 

0 

0 

.3 

4.000000 

5.886701 

3.321925 

2.000000 

.678076 

.113298 

0 

0 

0 

0 

0 

4.000000 

4.000000 

3.886711 

3.521923 

2.000000 

.678075 

.115298 

0 

0 

0 

0 

.5 

4.000000 

4.000000 

4.000000 

3.886701 

3.321923 

2.000000 

.678076 

.113298 

0 

0 

0 

.6 

4.000000 

4.000000 

4.000000 

4.000000 

3. 886711 

3.321923 

2.000000 

.678075 

.113298 

0 

0 

• 7 

4.000000 

4.000000 

4.000000 

4.000000 

4.000000 

3.886701 

3.321923 

2.000000 

.678076 

.113298 

0 

.8 

5.837920 

3.857920 

3.857920 

3.837920 

3.837920 

3.837900 

3.724633 

3.159844 

1.837926 

.515998 

0 

• 9 

2.991145 

2.991145 

2.991145 

2.991145 

2.991145 

2.991145 

2.991145 

2.877837 

2.315059 

1.018592 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.2 

-0.047438 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.1 

-.165594 

-.047438 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.165594 

-.047438 

0 

0 

0 

0 

0 

0 

0 

0 

.1 

.165594 

0 

.165594 

-.047438 

0 

0 

0 

0 

0 

0 

0 

.2 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

0 

0 

0 

0 

0 

• 3 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

0 

0 

0 

0 

.4 

0 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

0 

0 

0 

.5 

0 

0 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

0 

0 

.6 

0 

0 

0 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

0 

• 7 

0 

0 

0 

0 

0 

.047438 

.165594 

0 

-.165594 

-.047438 

0 

.8 

.679044 

.067904 

.067904 

.067904 

.067904 

.067904 

.115356 

.233502 

.067905 

-.097688 

0 

.9 

.249504 

.259504 

.249504 

.249504 

.249504 

.249504 

.249504 

.296940 

.414640 

.237672 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABLE  II.  - SECTION  LIFT  AND  MOMENT  COEFFICIENTS  Be,  AND  $c_ 

‘P  “P 

FOR  A UNIT  RATE  OF  ROLL 


^p  "'ifr  0cia  + rrr 

pCfflP  * "lTT  pc%  + ITT  PS0 


y h 

?cia(y/i) 

0Cjpo(y/l) 

(a)  31/e  - 5/7 

0 

5.714314 

0.306768 

0.109560 

-0.118325 

.1 

3.589100 

.047812 

.145258 

-.159784 

.2 

3.442694 

-.184467 

.180026 

-.201629 

.3 

3.273239 

-.387793 

.212207 

-.241915 

.4 

3.077662 

-.558969 

.240010 

-.278412 

.5 

2.851240 

-.695521 

.261463 

-.308527 

.6 

2.586448 

-.784948 

.273900 

-.328600 

.7 

2.270441 

-.823060 

.273661 

-.333866 

.8 

1.878086 

-.789604 

.254594 

-.315697 

• 9 

1.344769 

-.641916 

.202991 

-.255769 

1.0 

0 

0 

0 

0 

(b)  pl/c  - 1 

0 

4.000000 

0 

0 

0 

.1 

3.944613 

-.343532 

.025464 

-.025974 

.2 

3.837918 

-.631772 

.067905 

-.070621 

.3 

3.690902 

-.873882 

.116694 

-.123696 

.4 

3.503889 

-1.068886 

.166334 

-.179641 

.5 

3.273240 

-1.212206 

.212207 

-.233427 

.6 

2.991135 

-1.295676 

.249503 

-.279443 

• 7 

2.642985 

-1.305516 

.272287 

-.310407 

.8 

2.199256 

-1.216158 

.271624 

-.315084 

•9 

1.583274 

-.966580 

.229183 

-.270436 

1.0 

0 

0 

0 

0 

(c)  Bl/c  - 4/3 

0 

4.000000 

0 

0 

0 

.1 

4.000000 

-.400000 

0 

0 

.2 

4.000000 

-.800000 

0 

0 

.3 

3.970171 

-1.169877 

.014115 

-.014256 

.4 

3.837919 

-1.433309 

.067906 

-.069943 

• 5 

3.C33116 

-1,616485 

.153380 

-.140049 

.6 

3.355351 

-I.716783 

.197619 

-.211452 

.7 

2.991142 

-1.7195^3 

.249502 

-.271957 

.8 

2.507913 

-1.593429 

.275267 

-.305547 

.9 

1.817491 

-1.260633 

.250073 

-.282582 

1.0 

0 

0 

0 

0 
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TABLE  II.-  SECTION  LIFT  AND  MOMENT  COEFFICIENTS 


FOR  A UNIT  RATE  OF  ROLL  - Concluded 


%--rhSe‘a  + rhs%0 

%--rhsx  + r^%_ 


1 

iU 

y/i 

B«ta(y/i) 

Bclpo(y/l) 

sSo(y/0 

(d)  pi/c  - 

2 

4.000000 

0 

0 

0 

4.000000 

-.400000 

0 

0 

4.000000 

-.800000 

0 

0 

4.000000 

-1.200000 

0 

0 

4.000000 

-1.600000 

0 

0 

4.000000 

-2.000000 

0 

0 

3.837921 

-2.234846 

.067906 

-.069264 

3.503723 

-2.286270 

.166336 

-.172989 

2.991137 

-2.143407 

.249504 

-.264474 

2.199274 

-1.707722 

.271624 

-.293354 

0 

0 

0 

0 

(e)  pi/c  - 4 


4.000000 

0 

0 

0 

4.000000 

-.400000 

0 

0 

4.000000 

-.800000 

0 

0 

4.000000 

-1.200000 

0 

0 

4.000000 

-1.600000 

0 

0 

4.000000 

-2.000000 

0 

0 

4.000000 

-2.400000 

0 

0 

4.000000 

-2.800000 

0 

0 

3.837920 

-3.036384 

.067906 

-.068585 

2.991155 

-2.567274 

.249503 

-.256987 

0 

0 

0 

0 
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TABLE  III.-  SECTIONAL  LIFT  AND  MOMENT  COEFFICIENTS 


FOR  A UNIT  AILERON  DEFLECTION 


a £ 0. 14 

tgjl  - 1.0 

ca/c  “°*2 


y h 

3clg(y/z) 

^(y/z) 

(a)  pz/c  « 

5/7 

0 

-0.400000 

-0.160000 

.1 

-.646627 

-.251783 

.2 

-.760302 

-.301080 

.3 

-.800000 

-.320000 

.4 

-.800000 

-.320000 

.5 

-.800000 

-.320000 

.6 

-.800000 

-.320000 

.7 

-.800000 

-.320000 

.8 

-.742863 

-.292771 

• 9 

-.570248 

-.217640 

1.0 

0 

0 

(*>)  3l/c 

* 1 

0 

-0.400000 

-0.160000 

.1 

-.701015 

-.274665 

.2 

-.800000 

-.320000 

• 3 

-.800000 

-.320000 

.4 

-.800000 

-.320000 

.5 

-.800000 

-.320000 

.6 

-.800000 

-.320000 

.7 

-.800000 

-.320000 

.8 

-.800000 

-.320000 

• 9 

-.654649 

-.253370 

1.0 

0 

0 

y/i 

3e,6(y/l) 

PCngCy/z) 

(c)  pz/c  - 

4/3 

0 

-0.400000 

-0.160000 

.1 

-.749212 

-.295998 

.2 

-.800000 

-.320000 

.5 

-.800000 

-.320000 

.4 

-.800000 

-.320000 

.5 

-.800000 

-.320000 

.6 

-.800000 

-.320000 

.7 

-.800000 

-.320000 

.8 

-.800000 

-.320000 

• 9 

-.726623 

-.285314 

1.0 

0 

0 

(d)  pz/c 

£ 2 

0 

-0.400000 

-0.160000 

.1 

-.800000 

-.320000 

.2 

-.800000 

-.320000 

.3 

-.800000 

-.320000 

.4 

-.800000 

-.320000 

.5 

-.800000 

-.320000 

.6 

-.800000 

-.520000 

.7 

-.800000 

-.320000 

.8 

-.800000 

-.320000 

•9 

-.800000 

-.320000 

1.0 

0 

0 
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effectiveness, 

* 


.4 


0 .2  .4  .6  .8  1.0 

Pressure  rotio  , -jr 


(b)  M > 1.202. 


Figure  3*-  Variation  of  rolling  effectiveness  with  pressure  ratio  for 

constant  values  of  Much  number. 
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(c)  M * 1.338. 
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Figure  3--  Continued. 
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6.-  Variation  of  pressure  ratio  vlth  Mach  number  for  constant 
values  of  rolling  effectiveness. 
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Figure  7.-  Variation  of  rolling  effectiveness  with  Mach  number  at 

constant  altitude. 
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